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SUMMARY 


This report documents a FORTRAN program that provides a powerful and flexible 
tool for the linearization of aircraft models. The program LINEAR numerically 
determines a linear system model using nonlinear equations of motion and a user- 
supplied nonlinear aerodynamic model. The system model determined by LINEAR con- 
sists of matrices for both state and observation equations. The program has been 
designed to allow easy selection and definition of the state, control, and observa- 
tion variables to be used in a particular model. 


INTRODUCTION 


The program LINEAR was developed at the Dryden Flight Research Facility of 
NASA's Ames Research Center to provide a standard, documented, and verified tool 
to be used in deriving linear models for aircraft stability analysis and control 
law design. This development was undertaken to eliminate the need for aircraft- 
specific linearization programs common in the aerospace industry. Also, the lack 
of available documented linearization programs provided a strong motivation for 
the development of LINEAR? in fact, the only available documented linearization 
program that was found in an extensive literature search of the field is that of 
Kalviste (1980). 

Linear system models of aircraft dynamics and sensors are an essential part of 
both vehicle stability analysis and control law design. These models define the 
aircraft system in the neighborhood of an analysis point and are determined by the 
linearization of the nonlinear equations defining vehicle dynamics and sensors. 

This report describes LINEAR, a FORTRAN program that provides the user with a power- 
ful and flexible tool for the linearization of aircraft models. LINEAR is a program 
with well-defined and generalized interfaces to aerodynamic and engine models and is 
designed to address a wide range of problems without requiring program modification. 

The system model determined by LINEAR consists of matrices for both state and 
observation equations. The program has been designed to provide easy selection 
and definition of the state, control, and observation variables to be used in a 
particular model. Thus, the order of the system model is completely under user 
control. Further, the program provides the flexibility of allowing alternative 
formulations of both state and observation equations. 

LINEAR has several features that make it unique among the linearization programs 
common in the aerospace industry. The most significant of these features is flexi- 
bility. By generalizing the surface definitions and making no assumptions of sym- 
metric mass distributions, the program can be applied to any aircraft in any phase 
of flight except hover. The unique trimming capabilities, provided by means of a 
user-supplied subroutine, allow unlimited possibilities for trimming strategies and 
surface scheduling, which are particularly important for oblique-wing vehicles and 
aircraft having multiple surfaces affecting a single axis. The formulation of the 
equations of motion permit the inclusion of thrust-vectoring effects. The ability 
to select, without program modification, the state, control, and observation vari- 
ables for the linear models, combined with the large number of observation quan- 
tities availiabie, allows any analysis problem to be solved with ease. 


This report documents the use of the program LINEAR, defining the equations 
used and the methods employed to implement the program • The trimming capabilities 
of LINEAR are discussed from both a theoretical and an implementation perspective. 
The input and output files are described in detail. The user-supplied subroutines 
required for LINEAR are discussed, and sample subroutines are presented. 


NOMENCLATURE 


The units associated with the listed variables are expressed in a generalized 
system (given in parentheses). LINEAR will work equally well with any consistent 
set of units with two notable exceptions: the printed output and the atmospheric 

model. Both the printed output and the atmospheric model assume English units. 
Where applicable, quantities are defined with respect to the body axis system. 


Variables 


A 

A' 

a 

a n 

a n, i 
a nx 

a nx, i 

a ny 

a ny , i 

a nz 

a nz, i 

a x 

a y 


state matrix of the state equation, x = Ax + Bu; or, axial force (force) 

state matrix of the state equation, Cx = A'x + B*u 
speed of sound in air (length/sec) 
normal acceleration (g) 

normal acceleration, accelerometer not at center of gravity (g) 

x body axis accelerometer output, accelerometer at center of gravity (g) 

x body axis accelerometer output, accelerometer not at center of 
gravity (g) 

y body axis accelerometer output, accelerometer at center of gravity (g) 

y body axis accelerometer output, accelerometer not at center of 
gravity (g) 

z body axis accelerometer output, accelerometer at center of gravity (g) 

z body axis accelerometer output, accelerometer not at center of 
gravity (g) 

acceleration along the x body axis (g) 

acceleration along the y body axis (g) 

acceleration along the z body axis (g) 
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B control matrix of the state equation, x = Ax + Bu 

B* control matrix of the state equation, Cx = A'x + B'u 

b wingspan ( length ) 

C C matrix of the state equation, Cx = A’x + B'u; or, force or moment 

coefficient 

C D coefficient of drag 

C L coefficient of lift 

C£ coefficient of rolling moment 

C m coefficient of pitching moment 

C n coefficient of yawing moment 

C Y coefficient of sideforce 

CMj_ center of mass of ith engine 

c mean aerodynamic chord (length) 

D dynamic interaction matrix for state equation, x = Ax + Bu + Dv; or, 

drag force ( force ) 

D* dynamic interaction matrix for the state equation, Cx=A'x+B , u+D , v 

E dynamic interaction matrix for the observation equation, y = Hx + Fu + Ev 

E' dynamic interaction matrix for the observation equation, 

y=H'x+Gx+F'u+E , v 

E s specific energy (length) 

F feedforward matrix of the observation equation, y = Hx + Fu 

F total aerodynamic force acting at the aerodynamic center 

Fp engine thrust vector 

• 

F' feedforward matrix of the observation equation, y=H'x+Gx+F'u 
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fpa 

G 

g 

H 

H 1 
h 

he 

I 

Ie 

lx 

I xy 

X xz 

X Y 

*yz 

*z 

L 

£ 

M 

m 

m e 

N 

n 

p s 


flightpath acceleration (g) 

G matrix of the observation equation, y=H , x+Gx+F'u 
acceleration due to gravity (length/sec 2 ) 

observation matrix of the observation equation, y = Hx + Fu 

observation matrix of the observation equation, y=H , x+Gx+F'u 
altitude (length) 

angular momentum of engine rotor (mass-length 2 /sec ) 
aircraft inertia tensor (mass-length 2 ) 
rotational inertia of the engine (mass-length 2 ) 

x body axis moment of inertia (mass-length 2 ) 
x-y body axis product of inertia (mass-length 2 ) 
x-z body axis product of inertia (mass-length 2 ) 
y body axis moment of inertia (mass-length 2 ) 
y-z body axis product of inertia (mass-length 2 ) 
z body axis moment of inertia (mass-length 2 ) 

total body axis aerodynamic rolling moment (length-force); or, 
total aerodynamic lift (force) 

generalized length (length) 

Mach number; or, total body axis aerodynamic pitching moment 
( length-force) 

aircraft total mass (mass) 

mass of engine 

normal force (force); or, total body axis aerodynamic yawing moment 
( length-force) 

load factor 

specific power (length/sec) 

! 

i 

i 


4 


p 


Pa 


Pt 

q 

q 


<Ic 


Re, R$ 


Re 
Re • 
r 

S 

T 


T t 


u 

u 

V 

V c 

V e 

v 

V 

w 


w 

X 

Xip 

X 


roll rate (rad/sec); or, pressure (force/length 2 ) 
ambient pressure (force/length 2 ) 

total pressure (force/length 2 ) 
pitch rate (rad/sec) 

dynamic pressure (force/length 2 ) 
impact pressure (force/length 2 ) 
axis transformation matrices 
Reynolds number 

Reynolds number per unit length (length” 1 ) 
yaw rate (rad/sec) 

wing planform area (length 2 ) 

ambient temperature (K); or, total angular momentum (mass~length 2 /sec 2 ) ; 
or thrust (force) 

total temperature (K) 

velocity in x-axis direction (length/sec) 

control vector 

total velocity (length/sec) 

calibrated airspeed (knots) 

equivalent airspeed (knots) 

velocity in y-axis direction (length/sec) 

dynamic interaction vector 

vehicle weight (force) 

velocity in z-axis direction (length/sec) 
total force along the x body axis (force) 
thrust along the x body axis (force) 
state vector 
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Y 


sideforce (force) 


Y t thrust along the y body axis (force) 

y observation vector 

Z total force along the z body axis (force) 

Zrji thrust along the z body axis (force) 

a angle of attack (rad) 

3 angle of sideslip (rad) 

Y flightpath angle (rad) 

Ar displacement of aerodynamic reference point from center of gravity 

Ax displacement from center of gravity along x body axis (length) 

Ay displacement from center of gravity along y body axis (length) 

Az displacement from center of gravity along z body axis (length) 

6 a lateral trim parameter 

differential stabilator trim parameter 
6 e longitudinal trim parameter 

djfc Kronecker delta 

6 r directional trim parameter 

6 sb speed brake trim parameter 

thrust trim parameter 

6l incremental rolling moment (length-force) 

6m incremental pitching moment (length-force); or, incremental Mach 

<5n incremental yawing moment (length-force) 

6x incremental x body axis force (force) 

6y incremental y body axis force (force) 

6Z incremental z body axis force (force) 

e angle from the thrust axis of engine to the x-y body axis plane (rad) 
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C angle from the projection of F p onto the engine x-y plane to the local 

x axis (rad) 

X] angle from F p to the engine x-y plane (rad) 

0 pitch angle (rad) 

y coefficient of viscosity 

5 angle from the projection of F p onto the x-y body axis plane to the 

x body axis (rad) 

I p density of air (mass/length^ ) 

[ El total body axis rolling moment (length-force) 

Em total body axis pitching moment (length-force) 

En total body axis yawing moment (length-force) 

T torque from engines (length-force) 

Tg gyroscopic torque from engine (length-force) 

<j> roll angle (rad) 

<|>L tilt angle of acceleration normal to the flightpath from the vertical 

plane (rad) 

tJj heading angle (rad) 

u> total rotational velocity of the vehicle 

<D e engine angular velocity (rad/sec) 

Superscripts 

nondimensional version of variable 
. derivative with respect to time 

' T transpose of a vector or matrix 

Subscripts 

ar aerodynamic reference point 

D total drag 


E 


engine 


h altitude 

L total lift 

l rolling moment 

M Mach number 

m pitching moment 

max maximum 

min minimum 

n yawing moment 

o offset from center of gravity 

p roll rate 

q pitch rate 

r yaw rate 

s stability axis 

TP thrust point 

x along the x body axis 

Y sideforce 

y along the y body axis 

z along the z body axis 

0 standard day, sea level conditions; or, along the reference trajectory 

FORTRAN Variables 

AIX inertia about the x body axis 

AIXE engine inertia 

AIXY inertial coupling between the x and y body axes 

AIXZ inertial coupling between the x and z body axes 

AIY inertia about the y body axis 

AIYZ inertial coupling between the y and z body axes 

AIZ inertia about the z body axis 
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ALP 

ALPDOT 

AMCH 

AMSENG 

AMSS 

B 

■ BTA 
BTADOT 

( CBAR 
CD 

i 

CDA 

CODE 

CDO 

CDSB 

CL 

CLB 

CLDA 

CLDR 

CLDT 

CLFT 

CLFTA 

CLFTAD 

CLFTDE 

CLFTO 

CLFTQ 

CLFTSB 

CLP 


angle of attack 

time rate of change of angle of attack 
Mach number 

total rotor mass of the engine 

aircraft mass 

wingspan 

angle of sideslip 

time rate of change of angle of sideslip 
mean aerodynamic chord 
total coefficient of drag 

coefficient of drag due to angle of attack 
coefficient of drag due to symmetric elevator deflection 
drag coefficient at zero angle of attack 
coefficient of drag due to speed brake deflection 
total coefficient of rolling moment 

coefficient of rolling moment due to angle of sideslip 

coefficient of rolling moment due to aileron deflection 

coefficient of rolling moment due to rudder deflection 

coefficient of rolling moment due to differential elevator deflection 

total coefficient of lift 

coefficient of lift due to angle of attack 

coefficient of lift due to angle-of -attack rate 

coefficient of lift due to symmetric elevator deflection 

lift coefficient at zero angle of attack 

coefficient of lift due to pitch rate 

coefficient of lift due to speed brake deflection 

coefficient of rolling moment due to roll rate 
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CLR 

I 

CM 

I 

CMA 

CMAD 

CMDE 

i 

CMO 

CMQ 

CMSB 

CN 

CNB 

CNDA 

CNDR 

CNDT 

CNP 

CNR 

CY 

CYB 

i 

CYDA 

CYDR 

CYDT 

i 

DAS 

DC 

DELX 

DELY 

| DELZ 

DES 


coefficient of rolling moment due to yaw rate 
total coefficient of pitching moment 

coefficient of pitching moment due to angle of attack 

coefficient of pitching moment due to angle-of-attack rate 

coefficient of pitching moment due to symmetric elevator deflection 

pitching moment coefficient at zero angle of attack 

coefficient of pitching moment due to pitch rate 

coefficient of pitching moment due to speed brake deflection 

total coefficient of yawing moment 

coefficient of yawing moment due to sideslip 

coefficient of yawing moment due to aileron deflection 

coefficient of yawing moment due to rudder deflection 

coefficient of yawing moment due to differential elevator deflection 

coefficient of yawing moment due to roll rate 

coefficient of yawing moment due to yaw rate 

total coefficient of sideforce 

coefficient of sideforce due to sideslip 

coefficient of sideforce due to aileron deflection 

coefficient of sideforce due to rudder deflection 

coefficient of sideforce due to differential elevator deflection 

longitudinal trim parameter 

surface deflection and thrust control array 

displacement of the aerodynamic reference along the x body axis from 
the center of gravity 

displacement of the aerodynamic reference along the y body axis from 
the center of gravity 

displacement of the aerodynamic reference along the z body axis from 
the center of gravity 

lateral trim parameter 


10 


DRS 

DXTHRS 

EIX 

ENGOMG 

GMA 

H 

HDOT 

P 

PDOT 

PHI 

PHIDOT 

PSI 

i PSIDOT 

Q 

QBAR 

QDOT 

R 

RDOT 

S 

T 

| TDOT 

I 

1 THA 

( 

I 

THADOT 

THRSTX 

THRUST 

TLOCAT 


directional trim parameter 

distance between the center of gravity of the engine and the 
thrust point 

rotational inertia of the engine 
rotational velocity of the engine 
flightpath angle 
altitude 

time rate of change of altitude 
roll rate 

time rate of change of roll rate 
roll angle 

time rate of change of roll angle 
heading angle 

time rate of change of heading angle 
pitch rate 
dynamic pressure 

time rate of change of pitch rate 
yaw rate 

time rate of change of yaw rate 

wing area 

time 

time rate of change of time 
pitch angle 

time rate of change of pitch angle 

thrust trim parameter 

thrust generated by each engine 

location of the engine in the x, y , and z axes from the center of gravity 
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TVANXY orientation of the thrust vector in the x-y engine axis plane 

TVANXZ orientation of the thrust vector in the x-z engine axis plane 

UB velocity along the x body axis 

V velocity 

VB velocity along the y body axis 

VCAS calibrated airspeed 

VDOT time rate of change of total vehicle velocity 

VEAS equivalent airspeed 

WB velocity along the z body axis 

X position north from an arbitrary reference point 

XDOT time rate of change of north-south position 

XYANGL orientation of engine axis in x-y body axis plane 

XZANGL orientation of engine axis in x-z body axis plane 

Y position east from an arbitrary reference point 

YDOT time rate of change of east-west position 

PROGRAM OVERVIEW 


The program LINEAR numerically determines a linear system model using nonlinear 
equations of motion and a user-supplied nonlinear aerodynamic model, LINEAR is also 
capable of extracting linearized gross engine effects (such as net thrust, torque, 
and gyroscopic effects) and including these effects in the linear system model. The 
point at which this linear system model is defined is determined either by specifying 
the state and control variables or by selecting an analysis point on a trajectory, 
selecting a trim option, and allowing the program to determine the control vari- 
able and remaining state variables to satisify the trim option selected. 

Because the program is designed to satisfy the needs of a broad class of users, 
a wide variety of options has been provided. Perhaps the most important of these 
options are those that allow user specification of the state, control, and observa- 
tion variables to be included in the linear model derived by LINEAR. 

Within the program, the nonlinear equations of motion include 12 states repre- 
senting a rigid aircraft flying in a stationary atmosphere over a flat nonrotating 
earth. Thus, the state vector x is computed internally as 

x = [p q r V ot 3 <f> 0 h x y] T 
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The nonlinear equations used to determine the derivatives of the quantities are 
presented in the following section (Equations of Motion)* The internal control 
vector u can contain up to 30 controls* The internal observation vector y contains 
120 variables, including the state variables, the time derivatives of the state 
variables, the control variables, and a variety of other parameters of interest* 
Thus, within the program. 


y 

A 

II 

U T y? 

y T 

y 2 

yT yT yT yT yT yT^ T 

Y 3 y 4 y 5 y 6 y 7 y 8j 

yi = [ax 

a y a z a nx a ny 

a nz 

a n a nx,i a ny,i a nz,i a n,i 

y 2 = [a 

Re Re ' M 

q q c 

Pa 

T 

Pc /Pa Pt T T t V e V c ] 

U> 

II 

l—l 

• *' 

fpa y h 

h/57. 3 ; 

i T 


Y4 - [E s 

T 

Ps] 




y 5 = [L 

T 

D N A] 




Y6 = t u 

• • 

V w u V 

• T 

w] 




y? = [« #i e fi h fi h,i] T 

T 

Y8 = [T Ps qs r sl 


The equations defining these quantities are presented in the Observation 
Equations section* 


From the internal formulation of the 
state, control, and observation variables, 
the user must select the specific vari- 
ables desired in the output linear model 
(described in the Selection of State, Con- 
trol, and Observation Variables section)* 
Figure 1 illustrates the selection of var- 
iables in the state vector for a requested 
linear model. From the internal formula- 
tion on the right of the figure, the re- 
quested model is constructed, and the 
linear system matrices are selected in 
accordance with the user specification 
of the state, control, and observation 
variables. 


Output model Internal 

parameters parameters 



P 

q 

r 

V 

a 

P 

♦ 

e 

h 

x 


Figure 1. Selection of state vari- 
ables for linear model • 
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The linear model derived by LINEAR is determined at a specific analysis point. 
LINEAR allows this analysis point to be defined as a true steady-state condition 
on a specified trajectory (a point at which the rotational and translational accel- 
erations are zero) or as a totally arbitrary state on a trajectory. LINEAR provides 
the user with several options, described in detail in the Analysis Point Definition 
section. These analysis point definition options allow the user to trim the air- 
craft in wings-level flight, pushovers, pullups, level turns, or zero-sideslip ma- 
neuvers. Also included is a nontrimming option in which the user defines a totally 
arbitrary condition about which the linear model is to be derived. 

The linear system matrices are determined by numerical perturbation and are 
the first-order terms of a Taylor series expansion about the analysis point, as 
described in the Linear Models section. The formulation of the output system 
model is under user control. The user can select state equation matrices corre- 
sponding to either the standard formulation of the state equation, 

• 

x = Ax + Bu 

or the generalized equation, 

Cx = A'x + B'u 

The observation matrices can be selected from either of two formulations corre- 
sponding to the standard equation, 

y = Hx + Fu 

or the generalized equation. 


y=H , x+Gx+F , u 

In addition to the linear system matrices, LINEAR also computes the nondimen- 
sional stability and control derivatives at the analysis point. These derivatives 
are discussed in the Nondimensional Stability and Control Derivatives section. 

The input file for LINEAR is an ASCII file that defines the geometry and mass 
properties of the aircraft and selects various program options. Within this input 
file, the state, control, and observation vectors desired in the output linear model 
are defined, and the analysis point options are selected. The details of the input 
file are discussed in the Input Files section. 

The output of LINEAR is three files, one containing the linear system matrices 
and two documenting the options and analysis points selected by the user. The first 
is intended to be used with follow-on design and analysis programs. The other two 
contain all the information contained in the first file and also include the details 
of the analysis point and the nondimensional stability and control derivatives. 

These files are described in the Output Files section. 

To execute LINEAR, five user-supplied subroutines are required. These routines, 
discussed in the User-Supplied Subroutines section, define the nonlinear aerodynamic 
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model, the gross engine model, the gearing between the LINEAR trim inputs and the 
surfaces modeled in the aerodynamic model, and a model of the mass and geometry 
properties of the aircraft. The gearing model (fig. 2) defines how the LINEAR trim 
inputs will be connected to the surface models and allows schedules and nonstandard 
trimming schemes to be employed. This last feature is particularly important for 
oblique-wing aircraft. 


Inputs from 
/CTPARM/ 

DES 

DAS 

DRS 

THRSTX - 


UCNTRL 

(Gearing) 


Outputs to 
/ CONTRL/ 

► DC(1) 

► DC(2) 


DC(30) 


(Pilot stick, pedal, 
and throttle) 


(Surface deflections 
and power level setting) 


Figure 2 • Inputs to and outputs of the 
user-supplied subroutine UCNTRL • 


EQUATIONS OF MOTION 


The nonlinear equations of motion used in the linearization program are general 
six-degree-of-freedom equations representing the flight dynamics of a rigid aircraft 
flying in a stationary atmosphere over a flat nonrotating earth. The assumption of 
nonzero forward motion also is included in these equations; because of this assump- 
tion, these equations are invalid for vertical takeoff and landing or hover. These 
equations contain no assumptions of either symmetric mass distribution or aerody- 
namic properties and are applicable to asymmetric aircraft (such as oblique-wing 
aircraft) as well as to conventional symmetric aircraft. These equations of motion 
were derived by Etkin (1972), and the derivation will be detailed in a proposed NASA 
Reference Publication, "Derivation and Definition of a Linear Aircraft Model," by 
Eugene L. Duke, Robert F. Antoniewicz, and Keith D. Krambeer (in preparation). 

The following equations for rotational acceleration are used for analysis point 
definition: 


p = [(2 l)I 1 + (Zm)I 2 + (£n)I 3 - P 2 ( 1x2*2 - Ixyl3> 


+ pqdxzii 
+ q 2 (I yZ I 1 
- r 2 Uy Z Il 


I yz I 2 “ d z i 3 ^ “ P r ( I xy I 1 + D y I 2 “ 1 yz 1 3~> 
I xy I 3> _ g.r(D x Ii - Ixyl 2 + I xz I 3) 
I xz I 2)]/det I 
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q = [(£l)I 2 + (2M)l4 + (En)I 5 - P 2 (l xz l 4 - I xy I 5 ) 

+ P c t(I xz l2 — Iyz^4 - ~ P r ( ^xy I 2 + Py-^4 ~ Iyz^5^ 

+ q 2 (Iy Z l2 - I X yl5) ~ c l r (D x l2 - Ixy^4 + 

" r 2 (I yz I 2 " I xz l 4 )Vdet I 

r = [ (SL ) I 3 + (£M)I 5 + (SN)I 6 - P 2 (I xz l 5 - I X y I 6> 

+ pp(l xz l3 — Iyz^5 — '- ) z^6^ ~ P r ^xy^3 + ~ ^yz^6^ 

+ <J 2 (I y zl3 - I X yl6^ " c l r '(D X l3 “ ^xyl5 + ^xz^e) 

- r2(l yz i 3 - I xz l 5 )]/det I 

where 


det I 


- I x I y I z 


2 2 2 

2 I xy I xz I yz ~ ^x^z ~ I y I xz ” ^z^xy 


2 

*1 ~ x y*z “ lyz 

1 2 = ^xy^z + ^yz^xz 

1 3 = J xy^yz + *y*xz 

2 

1 4 - Ix X Z ~ *xz 


I 5 “ ^x^yz + I xy I xz 


x 6 ~ I x I y ” *xy 
D x - Iz - Iy 


Dy - I X - I z 


D z “ Iy - lx 


Here, the body axis rates are designated p, q, and r, corresponding to roll rate, 
pitch rate, and yaw rate. The total moments about the x, y, and z body axes (rol- 
ling, pitching, and yawing moments) are designated El, Em, and En, respectively. 
These total moments are the sums of all aerodynamic moments and powerplant-induced 
moments due to thrust asymmetries and gyroscopic torques. (The equations used 
to determine the change in moment coefficients due to the noncoincidence of the 
vehicle center of gravity and the reference point of the nonlinear aerodynamic model 
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are derived in appendix A. The equations defining the engine torque and gyroscopic 
contributions to the total moments are derived in appendix B. ) The body axis moments 
and products of inertia are designated I x , l y , l z , i xy/ i xz , and I yz . These moments 

and products of inertia are elements of the inertia tensor I , defined as 



~*xy 

-Ixzl 

>1 

X 

w 

l 

)! 

H 

I y 

1 

H 

N 

[f^XZ 

~*yz 

J 


To derive the state equation matrices for the generalized formulation, 

I Cx = Ax + Bu 

(where A and B are the state and control matrices of the state equation), the rota- 
tional accelerations are cast in a decoupled-axes formulation. The equations used 
to derive the linearized matrices are 



The translational acceleration equations used in the program LINEAR for both 
analysis point definition and perturbation are 
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V = [-D cos 3 + Y sin 0 + X T cos a cos 0 + Y T sin 0 + Z T sin a cos 0 

- mg (sin 0 cos a cos 0 - cos 0 sin <p sin 0 - cos 0 cos <)> sin a cos 0)]/m 

a = [-L + Z T cos a - X T sin a + mg(cos 0 cos cos a + sin 0 sin a)]/Vm cos 0 

+ q - tan 0 (p cos a + r sin a) 

0 = [D sin 0 + Y cos 0 - X T cos a sin 0 + Y T cos 0 - Z T sin a sin 0 

+ mg (sin 0 cos a sin 0 + cos 0 sin <f> cos 0 - cos 0 cos <f> sin a sin 0)]/Vm 

+ p sin a - r cos a 

where a, 0, 0, and <j> are angles of attack, sideslip, pitch, and roll, respectively; 
Xt, Yip, and Z T are thrust along the x, y, and z body axes; and D is drag force, 

g gravitational acceleration, L total aerodynamic lift, m total aircraft mass, 

V total velocity, and Y sideforce. 

The equations defining the vehicle attitude rates are 

4> = p + q sin (J) tan 0 + r cos 4> tan 0 

0 = q cos <p - r sin 4> 

4> = q sin <J) sec 0 + r cos 4) sec 0 

where 4 1 is heading angle. 

The equations defining the earth-relative velocities are 

h = V(cos 0 cos a sin 0 - sin 0 sin 4> cos 0 - cos 0 sin a cos <p cos 0) 

x = V[cos 0 cos a cos 0 cos ^ + sin 0 (sin 4> sin 0 cos ^ ~ cos 4> sin 40 

+ cos 0 sin a (cos 4> sin 0 cos ^ + sin 4 > s ^ n 

y = V[cos 0 cos a cos 0 sin \p + sin 0 (sin <j) sin 0 sin tp + cos <\> cos 4>) 

+ cos 0 sin a (cos 4> sin 0 sin 4^ - sin 4> cos 4^)1 

where h is altitude. 

OBSERVATION EQUATIONS 

The user-selectable observation variables computed in LINEAR represent a broad 
class of parameters useful for vehicle analysis and control design problems. These 
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variables include the state, time derivatives of state, and control variables. Also 
included are air data parameters, accelerations, flightpath terms, and other miscel- 
laneous parameters. The equations used to calculate those parameters are derived 
from a number of sources (Clancy, 1975; Dommasch et al., 1967; Etkin, 1972; Gainer 
and Hoffman, 1972; Gracy, 1980). Implicit in many of these observation equations is 
an atmospheric model. The model included in LINEAR is derived from the U.S. Standard 
Atmosphere ( 1 962 ) . 

The vehicle body axis accelerations constitute the set of observation variables 
that, except for state variables themselves, are most commonly used in the aircraft 
control analysis and design problem. These accelerations are measured in g units 
and are derived directly from the body axis forces defined in the previous section 
i for translational acceleration. The equations used in LINEAR for the body axis 
accelerations a x , ay, and a z are 

a x “ (X-ji - D cos ot + L sin a - gm sin 6)/g Q m 

ay = (Y t + Y + gm cos 0 sin <|))/gQm 

a z = (Zt - D sin a - L cos a + gm cos 0 cos <l>)/gQin 

where subscript 0 denotes standard day, sea level conditions. The equations for the 
outputs of the body axis accelerometers (denoted by subscript n) that are at vehicle 
center of gravity are 


a nx - ( X T - D cos a + L sin a)/g Q m 
a ny = (Yt + Y)/gQin 
a nz = ( Z T ~ D s ± n a ** L cos a)/gQm 
a n = (-Z T + D sin a + L cos ot)/g 0 m 

For orthogonal accelerometers that are aligned with the vehicle body axes but are 
not at vehicle center of gravity (denoted by subscript ,i), the following equations 
apply: 


a nx, i = a nx - t(q 2 + r 2 )x x - (pq - r)y x - (pr + q)z x ]/g 0 


*ny 


,i = a ny + [(pq + r)x y - (p 2 + r 2 )y y + (qr - p)z y ]/g ( 


a nz ,i = a nz + [(pr - q)x z + (qr + p)y z - (q 2 + p 2 )z z ]/g ( 
a n,i = a n - !<P r - q)*z + (qr + p)vz - (q 2 + p 2 > z z j/g 0 


where the subscripts x, y, and z refer to the x, y, and z body axes, respectively, 
and the symbols x, y, and z refer to the x, y, and z body axis locations of the 
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sensors relative to the vehicle center of gravity. Also included in the set of 
acceleration equations is load factor , n = L/W, where L is the total aerodynamic 
lift and W is the vehicle weight. 


The air data parameters having the greatest application to aircraft dynamics and 
control problems are the sensed parameters and the reference and scaling parameters. 
The sensed parameters are impact pressure q c , ambient or free-stream pressure p a , 

total pressure p t , ambient or free-stream temperature T, and total temperature T^. 

The selected reference and scaling parameters are Mach number M, dynamic pressure q, 
speed of sound a, Reynolds number Re, Reynolds number per unit length. Re 1 , and the 
Mach meter calibration ratio q c /p a . These quantities are defined as 



a 


Re 


pv£ 

U 


Re 


pv 

U 


q 


q 


c 



pv 2 

2 


[(1 .0 + 0.2M 2 ) 3 * 5 - 1 .o]p 

r,.»» ( - f 5 

|_ \5.6M 2 - 0.8/ 


0 


- 1.0 p 


/ (1.0 + 0.2M 2 ) 3 * 5 - 1.0 


. 2M 2 ( -- 5 -- 7 - 6 - M ~ Y’ 5 ]- 1.0 

\5.6M 2 - 0.8/ J 


(M < 1.0) 
(M > 1.0) 

(M < 1.0) 
(M > 1.0) 


P 


t 


p a + q 

U 


T t = T(1 .0 + 0.2M 2 ) 


where % is length, p pressure, p the density of the air, and y the coefficient of 
viscosity. Free-stream pressure, free-stream temperature, and the coefficient of 
viscosity are derived from the U.S. Standard Atmosphere (1962). 

Also included in the air data calculations are two velocities: equivalent 

airspeed V e and calibrated airspeed V c , both computed in knots. The calculations 

assume that internal units are in the English system. The equation used for 
equivalent airspeed is 
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V e - 17.17 Yq (lb/ft 2 ) 

which is derived from the definition of equivalent airspeed , 

V e * 



where p Q = 0.002378 slug/ft 3 and V e is converted from feet per second to knots. 
Calibrated airspeed is derived from the following definition of impact pressure: 


q 


c 


.2 


/l .0 + — 

Ve 2 V ' 5 - ll 


lA 7 * 

°P 0 ) J 


/ TT \ 9 



/M n 

5.76 

1 

( a n ) p 0 


9 

\ 0 / 

5.6 - 0.8 (a 0 /V c ) J 


2.5 


- Pn 


(V C < a Q ) 
(V c > a Q ) 


For the case where V c < a Q , the equation for V c is 


V c = 1479.116 



- 1 .0 


(V c < a Q ) 


Calibrated airspeed is found using an iterative process for the case where V c > a Q : 


V c 



,0 - 


1 .0 


7.0(V c /a 0 )‘ 


2.5 


(V C > a Q ) 


is executed until the change in V c from one iteration to the next is less than 
0.001 knots. 


Also included in the observation variables are the f lightpath-related parameters 
(described in app. D), including flightpath angle y, flightpath acceleration fpa. 


•• • 

vertical acceleration h, flightpath angle rate y, 


and (for lack of a better category 


in which to place it) scaled altitude rate h/57.3. The equations used to determine 
these quantities are 



h = a x sin 0 - a y sin $ cos 0 - a z cos <j> cos 0 
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•• • • 

V h - hV 

V Vv 2 - h 2 


Two energy-related terms are included in the observation variables: 
energy E s and specific power P s , defined as 

V 2 
2g 


specific 


E s “ 


h + 


P 


s 



VV 

+ — 

g 


The set of observation variables available in LINEAR also includes four force 
parameters: total aerodynamic lift L, total aerodynamic drag D, total aerodynamic 

normal force N , and total aerodynamic axial force A • These quantities are defined as 


L = qSC L 


D = qSC D 

N = L cos a + D sin a 
A = -L sin a + D cos a 

where C D and Cl are coefficients of drag and lift, respectively. 

Six body axis rates and accelerations are available as observation variables. 

These include the x body axis rate u, the y body axis rate v, and the z body axis 

• • • 

rate w. Also included are the time derivatives of these quantities, u, v, and w. 
The equations defining these quantities are 

u = V cos a cos 3 

v = V sin 3 

w = V sin a cos 3 

• /X T - gm sin 0 - D cos a + L sin a\ 

u = ^ ) + rv - qw 

• / y T + 9™ cos 9 sin 4> + Y\ 

V = l m / + pw - ru 

. / Zrp + gm cos 0 cos 4> — D sin a - L cos a\ 

w = ^ ) + qu - pv 

The final set of observation variables available in LINEAR is a miscellaneous 
collection of other parameters of interest in analysis and design problems. The 
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first group consists of measurements from sensors not located at the vehicle center 
of gravity. These represent angle of attack a angle of sideslip 3 altitude 

h f i , and altitude rate h^j_ measurements displaced from the center of gravity by some 


x, y, and z body axis 

distances. 

The equations used to compute these quantities are 

a i 

= a - 

/qx - 

py\ 

r 1 

\ v 

/ 


= e + 

/rx - 
\ V 

pz\ 


h^j_ = h + x sin 0 - y sin <j) cos 0 - z cos <\> cos 0 

• • • 

h ,i = h + 0(x cos 0 + y sin <(> sin 0 + z cos <j> sin 0) 


- 4> ( y cos <(> cos 0 - z sin <f> cos 0) 

The remaining miscellaneous parameters are total angular momentum T, stability 
axis roll rate p s/ stability axis pitch rate q S/ and stability axis yaw rate r s/ 
defined as 

T = T ” 2 I X yPq ” 2I xz P r + lyq 2 ” 2I yzq r + 

p s = p cos a + r sin a 

q s = q 

r s = -p sin a + r cos a 

SELECTION OF STATE , CONTROL, AND OBSERVATION VARIABLES 


The equations in the two preceding sections define the state and observation 
variables available in LINEAR. The control variables are defined by the user in the 
input file. Internally, the program uses a 12-state model, a 120-variable obser- 
vation vector, and a 30-parameter control vector. These variables can be selected 
to specify the formulation most suited for the specific application. The order and 
number of parameters in the output model is completely under user control. Figure 1 
illustrates the selection of variables for the state vector of the output model. 
However, it should be noted that no model order reduction is attempted. Elements of 
the matrices in the internal formulation are simply selected and reordered in the 
formulation specified by the user. 

The selection of specific state, control, and observation variables for the for- 
mulation of the output matrices is accomplished by alphanumeric descriptors in the 
input file. The use of these alphanumeric descriptors is described in the Input 
Files section. Appendix C lists the state variables and their alphanumeric descrip- 
tors. Appendix D lists the observation variables and their alphanumeric descriptors. 
The alphanumeric descriptors for the selection of control parameters to be included 
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in the observation vector are the control variable names defined by the user in the 
input file, as described in the Input Files section. 


LINEAR MODELS 


The linearized system matrices computed by LINEAR are the first-order terms of a 
Taylor series expansion about the analysis point (Dieudonne, 1978; Kwakernaak and 
Sivan, 1972; proposed NASA RP by Duke, Antoniewicz, and Krambeer, in preparation) 
and are assumed to result in a time-invariant linear system. The validity of this 
assumption is discussed in the Analysis Point Definition section. The technique 
employed to obtain these matrices numerically is a simple approximation to the par- 
tial derivative, that is, 


3f B f(X Q + Ax) - f(x 0 - Ax) 

3x 2 Ax 

where f is a general function of x, an arbitrary independent variable. The Ax 
may be set by the user but it defaults to 0.001 for all state and control parameters 
with the single exception of velocity V, where Ax is multiplied by a, the speed of 
sound, to obtain a reasonable perturbation size. 

From the generalized nonlinear state. 


• • 

TX = f(x, X, u) 

and observation equations, 

y = g( x, x, u) 

the program determines the linearized matrices for the generalized formulation of 
the system: 

C fix = A 1 fix + B' 5u 


where 


fiy = H 1 fix + G fix + F 1 fiu 


C 


A 1 


B' 


T - 


3f 

3x 

_3_f 

9u 


3f 

3x 


H 1 


3g 

3x 
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g -la 

3x 


F' 


= la 

3u 


with all derivatives evaluated along the nominal trajectory defined by the analysis 
point (x Q , x Q , u Q ) ; the state, time derivative of state, and control vectors can be 

expressed as small perturbations (<Sx, 6x, 6u) about the nominal trajectory, so that 

x = x Q + <5x 

• • • 

X = x Q + fix 

U S U Q + 511 

In addition to the matrices for this generalized system, the user has the option 
of requesting linearized matrices for a standard formulation of the system: 


where 


fix = A fix 4- B 6u 
fiy = H fix + F 6u 


A 


B 


T 


3f 


-1 


3f 


3x 


3x 


T 


8f 

3x 


-1 


3£ 

3u 


H 


la + la 

8x 3x 


T 



3f 

3x 


_ 3g 3g 
F = -r 3 - + — 

3u 3x 


3f 

3x 


-1 


3f 

3u 


with all derivatives evaluated along the nominal trajectory defined by the analysis 
point (x Q , x Q , u Q ). 

LINEAR also provides two nonstandard matrices, D and E, in the equations 


X = Ax + Bu + Dv 
y = Hx + Fu + Ev 


25 


or D* and E* in the equations 


Cx=A'x+B , u+D'v 


y 


=H'x+Gx+F'u+ 


E' v 


These dynamic interaction matrices include the effect of external forces and moments 
acting on the vehicle. The components of the dynamic interaction vector v are incre- 
mental body axis forces (6x, 6 y, 6z) and moments (6l, 6m, 6n): 


Thus , 


and 


v 


’fix" 

6y 

6z 

6l 

6m 

6n_ 


D 1 


3f_ 

3v 


E* 


3g 

3v 


D 


T 


3f 

3x 


-1 


3f 

3v 


E 



3g 

3x 


T 


3f ' 


-1 


3f 


3x 


3v 


These matrices allow the effects of unusual subsystems or control effectors to be 
easily included in the vehicle dynamics. 

The default output matrices for LINEAR are those for the standard system for- 
mulation, However, the user can select matrices for either generalized or standard 
state and observation equations in any combination. Internally, the matrices are 
computed for the generalized system formulation and then combined appropriately to 
accommodate the system formulation requested by the user. 


ANALYSIS POINT DEFINITION 


The point at which the nonlinear system equations are linearized is referred to 
as the analysis point. This can represent a true steady-state condition on the spec- 
ified trajectory (a point at which the rotational and translational accelerations 
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are zero; Perkins and Hage, 1949; Thelander, 1965) or a totally arbitrary state on a 
trajectory, LINEAR allows the user to select from a variety of analysis points. 
Within the program, these analysis points are referred to as trim conditions, and 
several options are available to the user. The arbitrary state and control option 
is designated NOTRIM, and in selecting this option the user must specify all nonzero 
state and control variables. For the equilibrium conditions, the user specifies a 
minimum number of parameters, and the program numerically determines required state 
and control variables to force the rotational and translational accelerations to 
zero. The analysis point options are described in detail in the following sections. 

For all the analysis point definition options, any state or control parameter 
may be input by the user. Those state variables not required to define the analysis 
point are used as initial estimates for the calculation of the state and control 
conditions that result in zero rotational and translational accelerations. As each 
state variable is read into LINEAR, the name is compared to the list of alternative 
state variables names listed in appendix C. All state variables except velocity 
I must be specified according to this list. Velocity can also be defined by specify- 
i ing Mach number (see alternative observation variable names in app. D). Appendix E 
lists analysis point definition identifiers that are recognized by LINEAR. 

\ 

It should be noted that the option of allowing the user to linearize the system 
equations about a nonequilibrium condition raises theoretical issues (beyond the 
scope of this report) of which the potential user should be aware. Although all the 
analysis point definition options provided in LINEAR have been found to be useful in 
the analysis of vehicle dynamics, not all the linear models derived about these anal- 
ysis points result in the time-invariant systems assumed in this report. However, 
the results of the linearization provided by LINEAR do give the appearance of being 
time invariant. 

The linearization process as defined in this report is always valid for some 
time interval beyond the point in the trajectory about which the linearization is 
done. However, for the resultant system to be truly time invariant, the vehicle 
must be in a sustainable, steady-state flight condition. This requirement is 
something more than merely a trim requirement, which is typically represented as 
• 

x(t) = 0, indicating that for trim, all the time derivatives of the state vari- 
ables must be zero. (This is not the case, however: Trim is achieved when the 

acceleration-like terms are identically zero; no constraints need to be placed on 

the velocity-like terms in x. Thus, for the model used in LINEAR, only p, q, r, V, 

a, and 5 must be zero to satisfy the trim condition.) The trim condition is 
achieved for the straight-and-level , pushover-pullup, level turn, thrust-stabilized 
turn, and beta trim options described in the following sections. In general, the 
no-trim and specific power analysis point definition options do not result in a trim 
condition. 

Of these analysis point options resulting in a trim condition, only the straight- 
and-level and level turn options force the model to represent sustainable flight 
conditions. In fact, only in the special case where the flightpath angle is zero 
does this occur for these options. 

► 

As previously stated, the linearization of a nonlinear model and its represen- 
tation as a time-invariant system are always valid for some time interval beyond the 
analysis point on the trajectory. This time interval is determined by several fac- 
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tors (such as trim and sustainable flight conditions) and ultimately by accuracy 
requirements placed on the representation . Thus, in using the linear models pro- 
vided by this program, the user should exercise some caution. 


Un trimmed 

For the untrimmed option, the user specifies all state and control variables 
that are to be set at some value other than zero. The number of state variables 
specified is entirely at the user’s discretion. If any of the control variables 
are to be nonzero, the user must specify the control parameter and its value. The 
untrimmed option allows the user to analyze the vehicle dynamics at any flight con- 
dition, including transitory conditions. 


Straight-and-Level Trim 

The straight-and-level trims available in LINEAR are in fact wings-level, 
constant-f lightpath-angle trims. Both options available for straight-and-level 
trim allow the user to specify either a flightpath angle or an altitude rate. 
However, since the default value for these terms is zero, the default for both 
types of straight-and-level trim is wings-level, horizontal flight. 

The two options available for straight-and-level trim require the user to 
specify altitude and either an angle of attack or a Mach number. If a specific 
angle of attack and altitude combination is desired, the user selects the "Mach- 
trim" option, which determines the velocity required for the requested trajectory. 
Likewise, the "alpha-trim" option allows the user to specify Mach number and alti- 
tude, and the trim routines in LINEAR determine the angle of attack needed for the 
r eques ted tra j ec tory . 

The trim condition for both straight-and-level options are determined within 
the following constraints: 


P = q = r = 0 
<f> = 0 

The trim surface positions, thrust, angle of sideslip, and either velocity or angle 
of attack are determined by numerically solving the nonlinear equations for the 
translational and rotational acceleration. Pitch attitude 0 is determined by itera- 
tive solution of the altitude rate equation. 

Pushover-Pullup 

The pushover-pullup analysis point definition options result in wings-level 
flight at n t 1 . For n > 1 , the analysis point is the minimum altitude point of a 

pullup when h = 0. For n < 1, this trim results in a pushover at the maximum alti- 
tude with h = 0. There are two options available for the pushover-pullup analysis 
point definition: (1) the "alpha- trim" option in which angle of attack is deter- 

mined from the specified altitude, Mach number, and load factor and (2) the "load- 
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factor-trim” option in which angle of attack, altitude, and Mach number are spec- 
ified and load factor is determined according to the constraint equations. 

The analysis point is determined at the specified conditions, subject to the 
following contraints: 

p = r = 0 

q = -q [mg(n - cos(0 - a)) - Z T cos a + Xrn sin a] 

Vm cos p 

<(> = 0 


• • 

The expression for q is derived from the a equation by setting a = 0 and (J> = 0; 0 is 

derived from the h equation. The trim surface positions, thrust, angle of sideslip, 
and either angle of attack or load factor are determined by numerically solving the 
nonlinear equations for the translational and rotational accelerations. 


Level Turn 


The level turn analysis point definition options result in non-wings-level , 
constant- turn-rate flight at n > 1 . The vehicle model is assumed to have sufficient 
excess thrust to trim at the condition specified. If thrust is not sufficient, trim 
will not result, and the analysis point thus defined will have a nonzero (in fact, 
negative) velocity rate. 

The level trim options available in LINEAR require the specification of an alti- 
tude and a Mach number. The user can then use either angle of attack or load factor 
to define the desired flight condition. These two options are referred to as "alpha- 
trim" and "load-factor-trim," respectively. For either option, the user may also 
request a specific flightpath angle or altitude rate. Thus, these analysis point 
definitions may result in ascending or descending spirals, although the default is 
for the constant-altitude turn. 


The constraint equations for the coordinated level turn analysis point defini- 
tions are derived by Chen (1981) and Chen and Jeske (1981). Using the requested 
load factor, the tilt angle of the acceleration normal to the flightpath from the 
vertical plane, <t>L, is calculated using the equation 


^L = ± tan 


-1 (n 2 “ cos 2 W /2 

cos y 


where the positive sign is used for a right turn and the negative sign is used for a 
left turn. From 4 >l, turn rate can be calculated as 


^ = — tan <)>l 
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Using these two definitions, the state variables can be determined: 


q = ^ sin 2 <() L 


-sin y sin 3 ± ( sin 2 y sin 2 3 - ^ 

sin 2 (f) L 


tan 4 > l cos 3 


^ sin y 

p s = - "cos 3 " q tan P 

p = p s cos a - r s sin a 

r = p s sin a + r s cos a 

9 = sin -1 ( -p/ip ) 

♦ ® tan"l(q/r) 


The trim surface positions, thrust, angle of sideslip, and either angle of attack or 
load factor are determined by numerically solving the nonlinear equations for the 
translational and rotational accelerations. 


Thrust-Stabilized Turn 

The thrust-stabilized turn analysis point definition results in a constant- 
throttle, non-wings-level turn with a nonzero altitude rate. The two options 
available in LINEAR are "alpha-trim” and " load-fac tor-trim . " These options allow 
the user to specify either the angle of attack or the load factor for the analysis 
point. The altitude and Mach number at the analysis point must be specified for 
both options. The user also must specify the value of the thrust trim parameter 
by assigning a value to the variable THRSTX in the input file after the trim has 
been selected. 

The constraint equations for the thrust-stabilized turn are the same as those 
for the level turn. However, for this analysis point definition, flightpath angle 
is determined by LINEAR. 


Beta Trim 

The beta trim analysis point definition results in non-wings-level, horizontal 

flight with heading rate $ = 0 at a user-specified Mach number, altitude, and angle 
of sideslip. This trim option is nominally at 1 g, but as 3 varies from zero, 
normal acceleration decreases and lateral acceleration increases. For an aero- 
dynamically symmetric aircraft, a trim to 3 = 0 using the beta trim option results 
in the same trimmed condition as the straight-and-level trim. However, for an aero- 
dynamically asymmetric aircraft, such as an oblique-wing vehicle, the two trim 
options are not equivalent. 
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The constraint equations used with the beta trim option are 

p = q = r = 0 

The trim surface positions, thrust, angle of attack, and bank angle are determined 
by numerically solving the nonlinear equations for translational and rotational 
acceleration. Pitch attitude 6 is derived from the equation for flightpath angle y 
with Y =0: 

q = tan ~l / sin 3 sin j> + cos 8 s in o t\ 

V cos 3 cos a / 


Specific Power 

The specific power analysis point definition results in a level turn at a user- 
specified Mach number, altitude, thrust trim parameter, and specific power. Unlike 
the other trim options provided in LINEAR, the specific power option does not, in 

general, attempt to achieve velocity rate V = 0. In fact, because the altitude rate 

• 

h = 0 and specific power is defined by 


P S 



the resultant velocity rate will be 


p s g 

V 


• • • • • 

However, the other acceleration-like terms (p, q, r, a, and 3) will be zero if the 
requested analysis point is achieved. 


The constraint equations used with the specific power analysis point definition 
option can be derived from the load factor tilt angle equation used in the level turn 
analysis point definition with y = 0: 

<t> L = ±tan _1 (n 2 - 1) 1y/2 


(where the positive sign is used for a right turn and the negative sign is used for 
a left turn). 


i i 

i 

i 

q 


tan <t>L 
sin <(>l cos 3 


s tan d> L cos3 
p g = -q tan 3 
p = p s cos a - r s sin a 


i 
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r = p s sin a + r s cos a 


8 = sin"" 1 \ 

\ V 

(j) = tan*” 1 — 
r 

The analysis point surface positions, load factor, angle of attack, and angle of 
sideslip are determined by numerically solving the nonlinear equations for transla- 
tional and rotational accelerations. 


NONDIMENS IONAL STABILITY AND CONTROL DERIVATIVES 


The nondimensional stability and control derivatives computed by LINEAR from the 
nonlinear aerodynamic model assume broadly formulated linear aerodynamic equations. 
These equations include effects of what are normally considered exclusively lateral- 
directional parameters in the longitudinal force and moment coefficient equations 
and, conversely, effects of longitudinal parameters in the lateral-directional equa- 
tions. The reason for this is two-fold: application to a larger class of vehicle 

types, such as oblique-wing aircraft, and computational ease. 

The nondimensional stability and control derivatives assume the following equa- 
tions for the aerodynamic force and moment coefficients: 


°t " C * 0 + V * V + C *h Sh * C *M S ” + 5i 


+ Cn P + C q + c r + C o + C 8 

*3 


c m ~ C m 


n 

i=1 1 


+ c mpP + + c m r ^ + C m^ a + 


n 


Cn = 


°n + c n a a + C np 3 + C nh Sh + C nM Sm + £ C n6 Si 
u i = 1 1 


+ CjipP + C nq q + C nr r + c n ‘a + C n £g 


n 


Cl - c L + c^a + C Lp e + C Lh Sh + C Lm 6 m + J C L6 . Si 
u i=1 x 
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C D = C D 0 + C D a a + C D3^ + C D h 6h + C D m 514 + Z 5 i 

A A 

+ C Dp P + C Dq q + C Dr r + C D -a + C D *g 

n 

C Y = Cv + C y a + Cyg3 + Cy h Sh + C Ym 6m + C Y$ . 5 i 

u i=1 1 

A A 

+ c Y p P + C Yq q + C Yr r + C Y *a + C Y ‘0 


where the stability and control derivatives have the usual meaning. 


l 


i 


c e 


x 


ac e 


with C e being an arbitrary force or moment coefficient (subscript l denoting rolling 
moment, m pitching moment, L total lift, D total drag, n yawing moment, and Y force 
along the y body axis) and x an arbitrary nondimensional variable (a denoting angle 
of attack, £ angle of sideslip, h altitude, M Mach number, p roll rate, q pitch rate, 
and r yaw rate). The rotational terms in the equations are nondimensional versions 
of the corresponding state variable with 


a 

P 


bp 

2V 


A 

q 


cq 

2V 


r 


br 

2V 


a 


ca 

2V 


3 


I 


b§_ 

2V 


| where b is wingspan and c is mean aerodynamic chord. The 6-^ in the summations are 

the n control variables defined by the user. The effects of altitude and velocity 
[ (Mach) are included in the derivatives with respect to those parameters and in the 
incremental multipliers 

6h = h - h Q 
6 m = M - M Q 
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where the subscript zero represents the current analysis point (x Q , x , u ), 
described in the Analysis Point Definition section. 

All stability derivatives are computed as nondimensional terms except the alti- 
tude and velocity parameters. The control derivatives are in whatever units are 
used in the nonlinear aerodynamic model. The derivatives with respect to velocity 
are multiplied by the speed of sound (at the analysis point altitude) to convert 
them to derivatives with respect to Mach number. Derivatives with respect to angle 
of attack and angle of sideslip can be obtained in units of reciprocal degrees. 

These derivatives are simply the corresponding nondimensional derivatives multiplied 
by 180.0/*Fr. 


INPUT FILES 


The LINEAR input file (defined in table 1) is separated into seven major sec- 
tions: case title and file selection information; project title; geometry and mass 

data for the aircraft; state, control, and observation variable definitions for the 
state-space model; trim parameter specification; additional control surfaces that 
may be specified for each case; and various test case specifications. All the input 


data can be written on one file or various files according to 
above and according to the input format defined in table 1 . 
is listed in appendix F. All the input records to LINEAR are 

TABLE 1 . — INPUT FORMAT FOR LINEAR 

the divisions specified 
An example input file 
written in ASCII form. 

Input record 

Format 

Title and file selection information 

Case title 

( 20A4 ) 

Input file names 

( 6A1 0 ) 

Project title 

Project title 

(20A4) 

Geometry and mass data 



s 

b c 

Weight 

(4F1 3.0) 

lx 

Iy I z 

*xz J xy ^yz 

(6F1 3.0) 

DELX 

DELY 

DELZ LOGCG 

(3F10.0, 1 2A4 ) 

a min 

a max 


(2F1 3.0) 
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TABLE 1. — Continued. 



Input record 

Format 


State, control, and observation variable definitions 

NUMSAT EQUAT STAB 

( I 1 0, A4, 1 IX, A4) 


STATE 1 

DRVINC 1 

( 5A4, FI 0. 0 ) 


STATE 2 

DRVINC 2 

( 5A4, FI 0. 0 ) 


STATE 3 

• 

• 

DRVINC 3 

• 

• 

( 5A4, FI 0. 0 ) 

• 

• 


• 

NUMSUR 

. 

• 

(110) 


CONTROL 1 

LOCCNT 1 CONVR 1 CNTINC 1 

( 5A4, 11 0, A4, 6X, 

F10.0) 

CONTROL 2 

LOCCNT 2 CONVR 2 CNTINC 2 

( 5A4, 11 0, A4, 6X, 

F10.0) 

CONTROL 3 

. 

. 

LOCCNT 3 CONVR 3 CNTINC 3 

• • * 

. . • 

( 5A4, 11 0, A4, 6X, 

. 

• 

FI 0.0) 

• 

NUMYVC 

. . . 

EQUAT 

• 

(110, A4) 


MEASUREMENT 

1 PARAM 1 (1 to 3) 

(5A4, 3F10.0) 


MEASUREMENT 

2 PARAM 2 (1 to 3) 

(5A4, 3F10.0) 


MEASUREMENT 

. 

• 

• 

3 PARAM 3 (1 to 3) 

. 

• 

. 

(5A4, 3F10.0) 

. 

• 

• 


Trim parameter specification 

<$ e <$ e <5, . 6, 6 r . 6 r 6 th . <5 th 

mm e max d nun d max r mm r max I * I1 min tn max 

(8F10.0) 


Additional surface specification 

NUMXSR 


(IX, 12) 


ADDITIONAL 

SURFACE 1 LOCCNT 1 CONVR1 

( 5A4, 11 0, A4 ) 


ADDITIONAL 

SURFACE 2 LOCCNT 2 CONVR 2 

{ 5A4, 11 0, A4 ) 
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TABLE 1. —Concluded. 


Input records 

Format 


Additional surface specification 

( continued ) 


ADDITIONAL SURFACE 3 LOCCNT 3 C0NVR3 

( 5A4, 11 0, A4 ) 



Test case specification 


ANALYSIS 

POINT 

DEFINITION 

OPTION 

(20A4) 

ANALYSIS 

POINT 

DEFINITION 

SUBOPTION 

(A4) 

VARIABLE 

1 

VALUE 1 


( 5A4, FI 5. 5 ) 

VARIABLE 

2 

VALUE 2 


( 5A4, FI 5. 5 ) 

VARIABLE 

• 

• 

3 

VALUE 3 

• 

• 


(5A4, FI 5.5) 

• 

• 

• 

NEXT 


• 


• 

( A4 ) 

ANALYSIS 

POINT 

DEFINITION 

OPTION 

(20A4) 

ANALYSIS 

POINT 

DEFINITION 

SUBOPTION 

(A4 ) 

VARIABLE 

1 

VALUE 1 


( 5A4, FI 5.5) 

VARIABLE 

2 

VALUE 2 


(5A4,F1 5.5) 

VARIABLE 

• 

* 

3 

VALUE 3 

• 

• 


( 5A4, FI 5. 5 ) 

• 

* 

• 

END 


• 


• 

(A4) 

There 

are seven input 

files for the batch linearizer. 

The first is a file con- 


taining the vehicle title and the names of the six data input files* The six data 
input files are 

1. project title, 

2. mass and geometry properties, 

3. state, control, and observation vectors, 
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4# trim parameter limits, 

5. trim definition cases, and 

6. additional surfaces to be set# 

This is the order in which they must be defined# However, the sixth file will be 
read before the fifth file. 


Case Title, File Selection Information, and Project Title 

There are two title records that need to be specified for LINEAR: a title for 
the individual test cases (case title) being analyzed and the name of the specific 
vehicle (project title)# These records are read with a 20A4 format and are sepa- 
rated by a file selection record# Both titles appear on each page of the line 
printer output file# The file selection record contains the names of the files 
from which the data are read. The data contained on the files specified by the six 
fields of the file selection record are shown in table 2. The input file names are 
written in character strings 10 columns long, and if not specified, the data are 
assumed to be on the same file as the first title record and the file selection 
record. The local name of the file containing these first two records must be 
attached to FORTRAN device unit 1. 


TABLE 2. — DEFINITION OF FILES SPECIFIED 
IN FILE SELECTION RECORD 


Field number 

Data contained on 
selected file 

1 

Project title 

2 

Geometry and mass data 

3 

State, control, and 
observation variable 
definitions 

4 

Trim parameter 
specifications 

5 

Test case specification 

6 

Additional surface 
definitions 


Geometry and Mass Data 

The geometry and mass data file consists of four records that can either fol- 
low the title and file selection records on unit 1 or be stored on a separate file 
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defined by the first of the input file name definitions on the file selection 
record. The geometry and mass data records define the geometry , mass properties, 
location of the aerodynamic reference with respect to the center of gravity, and 
angle-of-attack range for the vehicle model being analyzed. The first record of 
this set of five records is a vehicle title record in which the vehicle name and 
specific configuration can be defined to describe the information presented in the 
remaining records. 

The second and third records define the vehicle geometry, mass, and mass dis- 
tribution. The second record defines the wing planform area S, in units of length 
squared, the wingspan b, in units of length, the mean aerodynamic chord of the wing, 

c, in units of length, and the sea level weight of the vehicle, Weight, in units of 
mass-length per second squared. The third record defines the vehicle moments and 
products of inertia in units of mass-length squared. The fourth record defines the 
location of the vehicle center of gravity with respect to the aerodynamic reference 
point. This fourth record defines the offset of the aerodynamic reference point with 
respect to the center of gravity of the vehicle in the normal right-handed body axis 
reference system with the positive x axis forward and DELX, DELY, and DELZ in table 1 
representing the x, y, and z body axis displacements of the aerodynamic reference in 
units of length (see app. A). The fourth variable of the fourth record is an alpha- 
numeric variable, read using a 1 2A4 format, to specify if corrections due to an off- 
set of the aerodynamic reference with respect to the center of gravity are to be 
computed in LINEAR or in the user-supplied aerodynamic subroutine CCALC, described 
in the User-Supplied Subroutines section. The variable LOGCG defaults to a state 
that causes the aerodynamic reference offset calculations to be performed by sub- 
routines within LINEAR. Any of the following statements in the LOGCG field will 
cause LINEAR not to make these corrections: 

NO CG CORRECTIONS BY LINEAR 

CCALC WILL CALCULATE CG CORRECTIONS 

FORCE AND MOMENT CORRECTIONS CALCULATED IN CCALC 

The final record of this geometry and mass data set defines the angle-of-attack 
range for which the user-supplied nonlinear aerodynamic model (CCALC) is valid. 

These parameters, <* m ± n and ot max , specify the minimum and maximum values of angle of 

attack to be used for trimming the aircraft model. These parameters are in units of 
degrees. 


State, Control, and Observation Variable Definitions 

The state, control, and observation variables to be used in the output formula- 
tion of the linearized system are defined in records that either follow the last of 
the previously described sets of records on FORTRAN unit 1 or are stored on a sep- 
arate file defined by the second of the input file selection definitions on the file 
selection record. The number of records in the state, control, and observation 
variable definition set is determined by the number of such variables defined by 
the user. 

The states to be used in the output formulation of the linearized system are 
defined in the first set of records in the state, control, and observation variable 
definitions. The first record of this set defines the number of states to be used 
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(NUMSAT), the formulation of the output equation (EQUAT), and whether the nondimen- 
sional stability derivatives with respect to angle of attack and angle of sideslip 
are to be output in units of reciprocal radians or degrees. The variable EQUAT is 
read using an A4 format and is tested against the following list: 

NONSTANDARD 

NON-STANDARD 

GENERALIZED 

EXTENDED 


If EQUAT matches the first four characters of any of the listed words, the out- 
put formulation of the state equation is 


Cx = A ' x + B f u 

If EQUAT is read in as STANDARD or does not match the preceding list, then the 
default standard bilinear formulation of the state equation is assumed, and the 
output matrices are consistent with the equation 

x = Ax + Bu 

The variable STAB is also read using an A4 format and is compared with the following 
list: 


DEGREES 

DGR 

If STAB matches the first four characters of either of these words, the nondimen- 
sional stability derivatives with respect to angle of attack and angle of sideslip 
are printed in units of reciprocal degress on the printer file. Otherwise, these 
derivatives are printed in units of reciprocal radians. 

The remaining records of the state variable definition set are used to specify 
the state variables to be used in the output formulation of the linearized system 
and the increments to be used for the numerical perturbation described in the Linear 
Models section. The state variable names are checked for validity against the state 
variable alphanumeric descriptors listed in appendix C. If a name is not recognized, 
the variable is ignored and a warning message is written to the printer file. The 
increment to be used with any state variable (in calculating the A' and H' matrices) 
and the time derivative of that state variable (in calculating the C and G matrices) 
can be specified using the DRVINC variable. The units are the units of the state 
variable or its time derivative except for derivatives with respect to velocity. 

When DRVINC is specified for velocity, DRVINC specifies a Mach number increment. 

If DRVINC is not specified by the user, the default value of 0.001 is used. 

The next set of records in the state, control, and observation variable defini- 
tions are those defining the variables to be used in the control vector of the out- 
put model. The first record of this set defines the number of control parameters 
to be used (NUMSUR) . The remaining records define the names of these variables 
(CONTROL), their location (LOCCNT) in the common block /CONTRL/ (see the User- 
Supplied Subroutines section), the units associated with these control variables 
( CONVR ) , and the increments (CNTINC) to be used with these variables in determining 
the B' and F 1 matrices. 
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Because LINEAR has no default control variable names, the control variable names 
input by the user are used for subsequent identification of the control variables. 
Therefore, consistency in the use of control variable names is extremely important, 
particularly when the user attempts to establish control variable initial conditions 
when using the untrimmed analysis point definition option. 

The CONVR field in the control variable records is used to specify if the con- 
trol variables are given in degrees or radians. CONVR is read using an A4 format 
and is compared to the following list: 


DEGREES 

DGR 

RADIANS 

RAD 

If CONVR agrees with the first four characters of either of the first two names, it 
is assumed that the control variable is specified in units of degrees. If CONVR 
agrees with the first four characters of either of the last two listed names, it is 
assumed that the control variable is specified in units of radians. No units are 
assumed if CONVR does not agree with any of the listed names. When it is assumed 
that the control variable is specified in units of radians, the initial value of the 
control variable is converted to degrees before being written to the printer file. 

The variable CNTINC can be used to specify the increments used for a particular 
control surface when the B* and F 1 matrices are being calculated. It is assumed 
that the units of CNTINC agree with those used for the surface, and no unit conver- 
sion is attempted on these increments. If CNTINC is not specified for a particular 
surface, a default value of 0.001 is used. 

The final set of records in the state, control, and observation variable defini- 
tions pertain to the specification of parameters associated with the observation 
vector, observation equation, and observation parameters. The first record of this 
set defines the number of observation variables (NUMYVC) to be used in the output 
linear model and the formulation of the output equation (EQUAT). The remaining 
records in this set specify the variables to be included in the observation vector 
(MEASUREMENT) and any position information (PARAM) that may be required to compute 
the output model for a sensor not located at the vehicle center of gravity. 

The variable used to specify the formulation of the observation equation (EQUAT) 
is compared with the same list of names used to determine the formulation of the 
state equation. If it is determined that the generalized formulation is desired, 
the observation equation 


y=H'x+Gx+F'u 

is used. Otherwise, the standard formulation is assumed, and the form of the obser- 
vation equation used is 


y = Hx + Fu 

The records defining the observation variables to be used in the output for- 
mulation of the linear model contain a variable that includes the parameter name 
(MEASUREMENT) and three fields (PARAM) defining, when appropriate, the location of 
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the sensor relative to the vehicle center of gravity. The parameter name is com- 
pared with the list of observation variables given in appendix D. If the parameter 
name is recognized as a valid observation variable name, that observation variable 
is included in the formulation of the output observation vector. If the param- 
eter name is not recognized , an error message is printed and the parameter named 
is ignored. The three variables represented by PARAM(I), P ARAM ( 2 ) , and P ARAM ( 3 ) 
provide the x-axis , y-axis, and z-axis locations, respectively, of the measurement 
with respect to the vehicle center of gravity if the selected observation is one of 
the following: 


a x, i 

a y, i 

a z, i 

a n,i 


h 


i 


a 


i 


The unit associated with these variables is length. If the selected observation 
variable is not in the preceding list, the PARAM variables are not used. The sole 
exception to this occurs when Reynolds number is requested as an observation vari- 
able. In that case, PARAM ( 1 ) is used to specify the characteristic length. When 

no value is input for PARAM ( 1 ) , the mean aerodynamic chord c is used as the char- 
acteristic length. 


Trim Parameter Specification 

There is one record in the trim parameter specification set that is associated 
with the subroutine UCNTRL (described in the User-Supplied Subroutines section). 
This record specifies the limits to be used for the trim parameters <$ e , 6 a , S r , and 

representing the longitudinal, lateral, directional, and thrust trim parame- 
ters, respectively. The units associated with these parameters are defined by the 
implementation of UCNTRL. 


Additional Surface Specification 

The first record of this set of additional surface specifications defines the 
number of additional controls to be included in the list of recognized control names 
(NUMXSR). The purpose of defining these additional controls is to allow the user to 
set such variables as landing gear position, wing sweep, or flap position. Only the 
controls are defined in the additional surface specification records; actual control 
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positions are defined in the analysis point definition records • Because there may be 
no additional controls, these secondary trim parameter specification records may not 
be present. If such controls are defined, the records defining them will have the 
format specified in table 1. The control variable name (ADDITIONAL SURFACE), loca- 
tion (LOCCNT) in the common block /CONTRL/, and the units associated with this con- 
trol variable (CONVR) are specified for each additional control. 


Test Case Specification 

The test case specification records allow the user to define the flight condi- 
tion, or analysis point, at which a linear model is to be derived. Multiple cases 
can be included in the test case specification records. The final record in each 
set directs the program to proceed (NEXT) or to stop (END) execution. 

The first record of a test case specification set determines the analysis point, 
or trim, option to be used for the current case. The ANALYSIS POINT DEFINITION 
OPTION parameter is read in and checked against the list of analysis point defini- 
tion identifiers described in appendix E. The second record of a test case specifi- 
cation set, defining an analysis point definition suboption (ANALYSIS POINT DEFINI- 
TION SUBOPTION), will be read only if the requested analysis point definition option 
has a suboption associated with it. These suboptions are defined in the Analysis 
Point Definition section. The valid alphanumeric descriptors for these suboptions 
are described in appendix E. 

The remaining records in a test case specification set define test conditions or 
initial conditions for the trimming subroutines. These records consist of a field 
defining a parameter name (VARIABLE) and its initial condition (VALUE). These 
records may be in any order; however, if initial Mach number is to be defined, the 
initial altitude must be specified before Mach number if the correct initial veloc- 
ity is to be determined. The parameter names are checked against all name lists 
used within LINEAR. Any recognized state, time derivative of state, control, or 
observation variable will be accepted. Trim parameters also can be set in these 
records • 

In general, setting observation variables and time derivatives of the state 
variables has little effect. However, for some of the trim options defined in the 
Analysis Point Definition section, Mach number and load factor are used. The thrust 
trim parameter only affects the specific power trim. For the un trimmed option, the 
initial values of the state and control variables determine the analysis point com- 
pletely. For all other trim options, only certain states are not varied; all con- 
trols connected to the aerodynamic and engine model are varied. 


OUTPUT FILES 


There are three output files from LINEAR: a general-purpose analysis file, a 

printer file containing the calculated case conditions and the state-space matrices 
for each case, and a printer file containing the calculated case conditions only. 

The general-purpose analysis file contains the title for the cases being ana- 
lyzed; the state, control, and observation variables used to define the state-space 
model; and the state and observation matrices calculated in LINEAR. The C and G 


42 



matrices are printed only if the user has selected an appropriate formulation of 
the state and observation equations. The output for this file is assigned to 
FORTRAN device unit 15. An example of a general-purpose analysis file is pre- 
sented in appendix G f corresponding to the format of table 3. 


TABLE 3. — ANALYSIS FILE FORMAT 


Variable 

Format 

Title of the case 

{ 4A20 ) 

Title of the aircraft 

(4A20) 

Case number 

( //, 64X, 13 ) 

Number of states, controls, and outputs 

(17X,I2,22X,I2,22X,I3) 

State equation formulation 

(36X, 2A4) 

Observation equation formulation 

( 36X, 2A4 ) 
(////) 

State variable names, values, and units 

{IX, 5A4, 3X,E1 2.6, 2X,A20) 
(////) 

Control variable names, values, and units 

(IX, 5A4, 3X,E1 2.6, 2X,A20) 
(////) 

Dynamic interaction variable names and units 

(IX, 5A4, 1 7X, A20 ) 
(////) 

Output variable names, values, and units 

(1 X,5A4,3X,E1 2.6, 2X,A20) 

Matrix name 

( // » A8 , / ) 

A matrix 

( 5 ( E 1 3.6) ) 

Matrix name 

( // r A8 , / ) 

B matrix 

( 5 (El 3 . 6 ) ) 

Matrix name 

(//,A8,/) 

D matrix 

( 5 (El 3.6)) 

Matrix name 

( // , A8 , / ) 

C matrix (if general form chosen) 

( 5 ( E 1 3.6) ) 

Matrix name 

<//,A8,/) 

H matrix 

(5(E1 3.6) ) 
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TABLE 3. — Concluded 



Variable 

Format 

Matrix name 


(//,A8,/) 

F matrix 


(5(E1 3.6) ) 

Matrix name 


(//,A8,/) 

E matrix 


(5(E1 3.6) ) 

Matrix name 


(//,A8,/) 

G matrix (if 

general form chosen) 

( 5 (El 3.6) ) 


The titles are written on the first two records of the file in 80-character 
strings and are specified in LINEAR as the title of the vehicle and the title for 
the cases. The next record contains the number of the case as defined in LINEAR 
(999 cases maximum). 

The number of states, controls, and outputs used to define each case are written 
on the following record. The formulation of the state and observation equations 
are listed next, followed by the names and values of the states, controls, dynamic 
interaction variables, and outputs. These values are followed by the matrices that 
describe each case. 

The titles records appear only at the beginning of the file; all other records 
are duplicated for each subsequent case calculated in LINEAR. The matrices are 
written row-wise, five columns at a time, as illustrated in the following tabulation, 
which shows a system containing 7 states, 3 controls, and 11 outputs using the 
general state equation and standard observation equation. 

| 


Size of matrix Output formulation 


A 

- 

[ 7 

X 

7] 

A 

- 

[ 7 

X 

5] 








[ 7 

X 

2] 

B 

= 

[ 7 

x 

3] 

B 

= 

[ 7 

X 

3] 

D 


[ 7 

X 

6] 

D 

= 

[ 7 

X 

5] 








[ 7 

X 

1 ] 

C 

= 

[ 7 

X 

7] 

C 

= 

[ 7 

X 

5] 








[ 7 

X 

2] 

H 

= 

[11 

X 

7] 

H 

= 

[11 

X 

5] 








[11 

X 

2] 

F 

= 

[11 

X 

3] 

F 

= 

[11 

X 

3] 

E 

= 

[11 

X 

6] 

E 

= 

[11 

X 

5] 








[11 

X 

1 ] 
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The second output file, which is assigned to FORTRAN device unit 3, contains the 
values calculated in LINEAR describing each case. The first section of this file 
contains a listing of the input data defining the aircraft's geometry and mass prop- 
erties, variable names defining the state-space model, and various control surface 
limits characteristic of the given aircraft. Appendix H presents an example of this 
output file. 

The second section of this file contains the trim conditions of the aircraft at 
the point of interest. These conditions include the type of trim being attempted, 
whether trim was achieved, and parameters defining the trim condition. 

The values for the variables of the state-space model at the trim condition are 

• • • • • • 

also printed. If trim was not achieved, p, q, r, V, 3, and a (calculated from the 
equations of motion) and the force and moment coefficients are printed. Changes 
in the geometry and mass properties are also printed. 

The third section of this output file contains the nondimensional stability and 
control derivatives for the trim condition calculated. The static margin of the 
aircraft at the given flight condition is also printed. 

The final section of this output file contains the state and observation matri- 
ces for the given flight condition. The formulation of the state equations and only 
| the terms of the matrices chosen by the user to define the model are printed. A 
subset of this output file containing only the trim conditions is assigned to 
FORTRAN device unit 2. 

The third output file, which is assigned to FORTRAN device unit 2, contains the 
trim conditions of the aircraft at the point of interest. These conditions include 
I the type of trim being attempted, whether trim was achieved, parameters defining the 
I trim condition, and the static margin of the aircraft at the given flight condition. 
Appendix H presents an example of this file. 


USER-SUPPLIED SUBROUTINES 


There are five subroutines that must be supplied by the user to interface LINEAR 
with a specific aircraft's subsystem models: ADATIN, CCALC, IFENGN, UCNTRL , and 

MASGEO. The first two subroutines constitute the aerodynamic model. The subroutine 
IFENGN is used to provide an interface between LINEAR and any engine modeling rou- 
tines the user may wish to incorporate. UCNTRL converts the trim parameters used by 
LINEAR into control surface deflections for the aerodynamic modeling routines. The 
subroutine MASGEO allows the user to define the mass and geometry properties of the 
vehicle as a function of flight condition or control setting. The use of these sub- 
routines is illustrated in figure 3, which shows the program flow and the interaction 
of LINEAR with the user-supplied subroutines. These subroutines are described in 
detail in the following sections. Examples of these subroutines are given in appen- 
dix I. 
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Main program 


User-supplied 

subroutines 



Figure 3 • Program flow diagram showing communication 
with user-supplied subroutines . 


Aerodynamic Model 

It is assumed that the aerodynamic models consist of two main subroutines, 
ADATIN and CCALC. ADATIN is used to input the basic aerodynamic data from remote 
storage and can also be used to define aerodynamic data. CCALC is the subroutine 
that uses these aerodynamic data, the state variables, and the surface positions 
to determine the aerodynamic coefficients. Either routine may call other subrou- 
tines to perform related or required functions; however, from the point of view 
of the interface to LINEAR, only these two subroutines are required for an aero- 
dynamic model. 
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Externally, ADATIN has no interface to the program LINEAR* The subroutine 
is called only once when the aerodynamic data are input or defined. The calling 
program must provide ADATIN with the input devices it requires, but no other accom- 
modation is necessary. The aerodynamic data are communicated from ADATIN to CCALC 
through named common blocks that occur in only these two routines. 

The interface between CCALC and the calling program is somewhat more compli- 
cated. However, the interface is standard, and this feature provides a framework 
about which a general-purpose tool can be built. This interface consists of several 
named common blocks that are used to pass state variables, air data parameters, sur- 
face positions, and force and moment coefficients between CCALC and the calling 
program. CCALC is executed whenever new aerodynamic coefficients are required (for 
example, once per frame for a real-time simulation). 

The main transfer of data into the subroutine CCALC is through six named common 
blocks. These common blocks contain the state variables, air data parameters, and 
surface positions. The transfer of data from CCALC is through a named common block 
containing the aerodynamic force and moment coefficients. The details of these com- 
mon blocks follow. 

The common block /DRVOUT/ contains the state variables and their derivatives with 
respect to time. The structure of this common block is as follows: 

COMMON /DRVOUT/ T , 


P 

, Q 

, R 

t 

V 

, ALP 

, BTA 

r 

THA 

, PSI 

, PHI 

r 

H 

, x 

, Y 

r 


TDOT , 

PDOT , QDOT , RDOT , 

VDOT , ALPDOT, BTADOT , 

THADOT, PSIDOT, PHIDOT, 

HDOT , XDOT , YDOT 

The body axis rates p, q, and r appear as P, Q, and R, respectively. Total veloc- 
ity is represented by the variable V, altitude by H; angle of attack (ALP), angle 
of sideslip ( BTA ) , and their derivatives with respect to time (ALPDOT and BTADOT, 
respectively) are also contained within this common block. 

The common block /SIMOUT/ contains the main air data parameters required for 
the function generation subroutine. The variables in this common block are 
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COMMON /SIMOUT/ AMCH, QBAR, GMA, DEL, 


UB , VB , WB , 

VEAS , VCAS 

Mach number and dynamic pressure are the first entries in the common block, sym- 
bolized by AMCH and QBAR, respectively. The body axis velocities u, v, and w are 
included as UB, VB, and WB, respectively. 

The /CONTRL/ common block contains the surface position information. The exact 
definition of each of the elements used for a particular aerodynamic model is deter- 
mined by the implementer of that model. The structure of the common block /CONTRL/ 
is as follows: 

COMMON /CONTRL/ DC (30) 

The common block /DATAIN/ contains geometry and mass information. 


COMMON /DATAIN/ S , 

B 

CBAR, 

AMSS 

AIX , 

AIY , 

AIZ , 

AIXZ 

AIXY , 

AIYZ, 

AIXE 



The first three variables in the common block, S, B, and CBAR, represent wing area, 
wingspan, and mean aerodynamic chord, respectively. The vehicle mass is represented 
by AMSS. 

The information on the displacement of the aerodynamic reference point of the 
aerodynamic data with respect to the aircraft center of gravity is contained in the 
/CGSHFT/ common block: 


COMMON /CGSHFT/ DELX, DELY, DELZ 

The displacements are defined along the vehicle body axis with DELX, DELY, and DELZ 
representing the displacements in the x, y, and z axes, respectively. 

The common block /SIMACC/ contains the accelerations, accelerometer outputs, and 
normal accelerometer output at the center of gravity of the aircraft. 

COMMON /SIMACC/ AX , AY , AZ , 

ANX, ANY, ANZ, 

AN 

The output common block /CLCOUT/ contains the variables representing the aerody- 
namic moment and force coefficients: 

COMMON /CLCOUT/ CL, CM, CN, CD, CLFT, CY 
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The variables CL, CM, and CN are the symbols for the coefficients of rolling moment 
C£, pitching moment C m , and yawing moment C n , respectively; these terms are body 
axis coefficients. The stability axis forces are represented by CD (coefficient 
of drag C D ), CLFT (coefficient of lift Cl), and CY (sideforce coefficient Cy). 


Control Model 

The program LINEAR attempts to trim the given condition using control inputs 
similar to those of a pilot: longitudinal stick, lateral stick, rudder, and throt- 

tle. The UCNTRL subroutine utilizes these trim output control values to calculate 
actual surface deflections and power level angles for the given aircraft (fig. 2). 
The location of each surface and power level angle in the /CONTRL/ common block is 
specified by the user in the input file (maximum of 30 surfaces). The limits for 
the control parameters in pitch, roll, yaw, and thrust are user specified (see Trim 
Parameter Specification in the Input Files section) and must agree in units with the 
calculations in CCALC. 

The common block /CTPARM/ contains the four trim parameters that must be related 
to surface deflections in the subroutine UCNTRL: 

COMMON /CTPARM/ DES, DAS, DRS , THRSTX 

The output from UCNTRL is through the common block /CONTRL/ described previously in 
the Aerodynamic Model section. The variables DES, DAS, DRS, and THRSTX are used to 
trim the longitudinal, lateral, directional, and thrust axes, respectively. 

For an aircraft using feedback (such as a statically unstable vehicle), state 
variables and their derivatives are available in the common block /DRVOUT/. If 
parameters other than state variables and their time derivatives are required for 
feedback, the user may access them using the common block /OBSERV/ described in the 
Mass and Geometry Model section of this report. 


Engine Model 

The subroutine IFENGN computes individual engine parameters used by LINEAR to 
calculate force, torque, and gyroscopic effects due to engine offset from the cen- 
terline. The parameters for each engine (maximum of four engines) are passed 
through common /ENGSTF/ as follows: 

COMMON /ENGSTF/ THRUST(4), TLOCAT(4,3), XYANGL ( 4 ) , XZANGL(4), 

TVANXY ( 4 ) , TVANXZ (4), DXTHRS ( 4 ) , 

EIX (4), AMSENG (4), ENGOMG ( 4 ) 

The variables in this common block correspond to thrust (THRUST); the x, y, and z 
body axis coordinates of the point at which thrust acts (TLOCAT) ; the orientation of 
the thrust vector in the x-y body axis plane (XYANGL), in degrees; the orientation 
of the thrust vector in the x-z body axis plane ( XZANGL ) , in degrees; the orienta- 
tion of the thrust vector in the x-y engine axis plane (TVANXY), in degrees; the 
orientation of the thrust vector in the x-z engine axis plane (TVANXZ), in degrees; 
the distance between the center-of-gravity of the engine and the thrust point 
(DXTHRS), measured positive in the negative x engine axis; the rotational inertia of 
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the engine (EIX); mass (AMSENG); and the rotational velocity of the engine (ENGOMG ) • 
The variables are all dimensioned to accommodate up to four engines. 

While the common block structure within LINEAR is designed for engines that do 
not interact with the vehicle aerodynamics, propeller-driven aircraft can be easily 
modeled by communicating the appropriate parameters from the engine model (IFENGN) 
to the aerodynamic model (CCALC) . 


Mass and Geometry Model 

The subroutine MASGEO allows the user to vary the center-of-gravity position and 
vehicle inertias automatically. Nominally, this subroutine must exist as one of the 
user subroutines, but it may be nothing more than RETURN and END statements. MASGEO 
is primarily for variable-geometry aircraft (such as an oblique-wing or variable- 
sweep configurations) or for modeling aircraft that undergo significant mass or iner- 
tia changes over their operating range. The center-of-gravity position and inertias 
may be functions of flight condition or any surface defined in the /CONTRL/ common 
block. Changes in center-of-gravity position are passed in the /CGSHFT/ common 
block, and inertia changes are passed in the /DATAIN/ common block. 

Care must be taken when using the subroutine MASGEO in combination with 
selecting an observation vector that contains measurements away from the center of 
gravity. If measurements are desired at a fixed location on the vehicle, such as a 
sensor platform or nose boom, the moment arm to the new center-of-gravity location 
must be recomputed as a result of the center-of-gravity shift for accurate results. 
This can be accomplished by implementing the moment arm calculations in one of the 
user subroutines and passing the new moment arm values through the /OBSERV/ common 
block: 


COMMON /OBSERV/ OBVEC(120), PARAM(120,6) 

The common block /OBSERV/ allows the user to access all the observation vari- 
ables during trim as well as to pass parameters associated with the observations 
back to LINEAR. The common block /OBSERV/ contains two single-precision vectors: 
OBVEC (120), and PARAM( 120,6). A list of the available observations and parameters 
is given in table 4. Access to the observation variables allows the user to imple- 
ment trim strategies that are functions of observations, such as gain schedules and 
surface management schemes. The parameters associated with the observations are 
used primarily to define the moment arm from the center of gravity to the point at 
which the observation is to be made. If the user subroutine MASGEO is used to 
change the center-of-gravity location and observations are desired at fixed loca- 
tions on the vehicle, then the moment arm from the new center-of-gravity location to 
the fixed sensor location ((x, y, z), in feet) must be computed in one of the user 
subroutines and passed back in the first three elements of the PARAM vector asso- 
ciated with the desired observation (PARAM(n = 1 to 3), where n is the number of 
the desired observation) . Additional information on observations and parameters 
can be found in the State, Control, and Observation Variable Definitions section. 
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TABLE 4. — OBSERVATION VARIABLES AVAILABLE IN THE USER-SUPPLIED SUBROUTINES 

USING THE OBVEC ARRAY 


Location (index) 
in OBVEC 

Variable 


State Variables 

1 

Roll rate 


2 

Pitch rate 


3 

Yaw rate 


4 

Velocity 


5 

Angle of attack 


6 

Angle of sideslip 


7 

Pitch attitude 


8 

Heading 


9 

Roll attitude 


10 

Altitude 


1 1 

Displacement north 


1 2 

Displacement east 


Derivatives of state variables 

1 3 

Roll acceleration 


14 

Pitch acceleration 


15 

Yaw acceleration 


16 

Velocity rate 


17 

Angle-of-attack rate 


1 8 

Angle-of-sideslip rate 


19 

Pitch attitude rate 


20 

Heading rate 


21 

Roll attitude rate 


22 

Altitude rate 


23 

Velocity north 


24 

Velocity east 


Accelerations 

25 

x body axis acceleration 


26 

y body axis acceleration 


27 

z body axis acceleration 


28 

x body axis accelerometer 

at vehicle 


center of gravity 


29 

y body axis accelerometer 

at vehicle 


center of gravity 


30 

z body axis accelerometer 

at vehicle 


center of gravity 


31 

Normal acceleration 


32 

x body axis accelerometer 

not at vehicle 


center of gravity 


33 

y body axis accelerometer 

not at vehicle 


center of gravity 
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TABLE 4. — * Continued. 


Location (index) 
in OBVEC 

Variable 


Accelerations (continued) 

34 

z body axis accelerometer not at vehicle 
center of gravity 

99 

Normal accelerometer not at vehicle 
center of gravity 

98 

Load factor 


Air data parameters 


91 

Speed of sound 

37 

Reynolds number 

1 03 

Reynolds number per unit length 

35 

Mach number 

36 

Dynamic pressure 

56 

Impact pressure 

55 

Ambient pressure 

57 

Impact/ambient pressure ratio 

58 

Total pressure 

59 

Temperature 

60 

Total temperature 

92 

Equivalent airspeed 

93 

Calibrated airspeed 

Flightpath-related parameters 

39 

Flightpath angle 

38 

Flightpath acceleration 

40 

Flightpath angle rate 

43 

Scaled altitude rate 

Energy-related terms 

46 

Specific power 

47 

Specific energy 


Force parameters 


94 

Lift force 

95 

Drag force 

96 

Normal force 

97 

Axial force 


Body axis parameters 


52 

x body 

axis 

velocity 

53 

y body 

axis 

velocity 

54 

z body 

axis 

velocity 
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TABLE 4. — Concluded. 


Location (index) 
in OBVEC 

Variable 


Body axis parameters (continued) 

100 

Rate of change of velocity 
in x body axis 

101 

Rate of change of velocity 
in y body axis 

102 

Rate of change of velocity 
in z body axis 

Miscellaneous 

measurements not at vehicle center of gravity 

44 

Angle of attack not at vehicle center 
of gravity 

45 

Angle of sideslip not at vehicle center 
of gravity 

41 

Altitude instrument not at vehicle 
center of gravity 

42 

Altitude rate instrument not at vehicle 
center of gravity 


Other miscellaneous parameters 


48 

49 

50 

51 

Vehicle total angular momentum 
Stability axis roll rate 
Stability axis pitch rate 
Stability axis yaw rate 

Control surface parameters 

61 to 90 

Control surfaces DC(1) to DC(30) 

Trim parameters 

1 04 

Longitudinal trim parameter 

105 

Lateral trim parameter 

106 

Directional trim parameter 

1 07 

Thrust trim parameter 


CONCLUDING REMARKS 


The FORTRAN program LINEAR was developed to provide a flexible, powerful, and 
documented tool to derive linear models for aircraft stability analysis and control 
law design. This report discusses LINEAR from the perspective of a potential user, 
defining the nonlinear equations from which the linear model is derived and describ- 
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ing the interfaces to user-supplied subroutines and input files. The output from 
LINEAR is also described. Examples of the user-supplied subroutines are presented 
in the appendixes. A microfiche listing of the program for a VAX 11/750 with the 
VMS operating system is included with this report. (Revisions to the microfiche 
supplement are given in app. J. ) 

The program LINEAR numerically determines a linear system model using nonlinear 
equations of motion and a user-supplied nonlinear aerodynamic model. LINEAR is 
also capable of extracting linearized engine effects (such as net thrust, torque, 
and gyroscopic effects) and including these effects in the linear system model. 

The point as which this linear system model is defined is determined either by 
completely specifying the state and control variables or by specifying an analysis 
point on a trajectory, selecting a trim option, and directing the program to deter- 
mine the control variables and remaining state variables. 

The system model determined by LINEAR consists of matrices for both state and 
observation equations. The program has been designed to provide an easy selection 
and definition of the state, control, and observation variables to be used in a 
particular model. Thus, the order of the system model is completely under user 
control. Further, the program provides the flexibility of allowing alternative 
formulations of both state and observation equations. 

LINEAR has several features that make it unique among the linearization programs 
common in the aerospace industry. The most significant of these features is flexi- 
bility. By generalizing the surface definitions and making no assumptions of sym- 
metric mass distributions, the program can be applied to any aircraft in any phase 
of flight except hover. The unique trimming capability, implemented through a user- 
supplied subroutine, allows unlimited possibilities of trimming strategies and sur- 
face scheduling, which are particularly important for oblique-wing vehicles and 
aircraft having multiple surfaces affecting a single axis. The formulation of the 
equations of motion permit the inclusion of thrust- vector ing effects. The ability to 
select, without program modification, the state, control, and observation variables 
for the linear models, combined with the large number of observation quantities 
available, allows any analysis problem to be solved with ease. 

This report documents the use of the program LINEAR, defining the equations 
used and the methods employed to implement the program. The trimming capabilities 
of LINEAR are discussed from both theoretical and implementation perspectives. 

The input and output files are described in detail. The user-supplied subroutines 
required for LINEAR are discussed, and sample subroutines are presented. 


National Aeronautics and Space Administration 
Ames Research Center 
Dryden Flight Research Facility 
Edwards, California, March 6, 1985 
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APPENDIX A: CORRECTION TO AERODYNAMIC COEFFICIENTS FOR A CENTER OF GRAVITY 

NOT AT THE AERODYNAMIC REFERENCE POINT 


The point on the vehicle at which the nonlinear force and moment coefficients 
are defined is referred to as the aerodynamic reference point. All aerodynamic 
effects are modeled at this aerodynamic reference point. Thus, when this point 
and the vehicle center of gravity are not coincident, the forces acting at the 
aerodynamic reference point effectively induce moments that act incrementally on 
the moments defined at the aerodynamic reference point by the nonlinear aerody- 
namic model. 

The total aerodynamic moment M acting at the vehicle center of gravity is 
defined as 


M = M ar -H Ar x F 


, where 

M ar = f L ar M ar N arl T 

is the total aerodynamic moment acting at the aerodynamic reference point (denoted 
by subscript ar) of the vehicle, 

Ar = [Ax Ay Az] T 

is the displacement of the aerodynamic reference point from the vehicle center of 
j gravity, and 

I F = [X Y Z]T 

: is the total aerodynamic force acting at the aerodynamic center, where X, Y, and Z 

\ are total forces along the x, y (sideforce), and z body axes. Thus, 

L ar + Ay Z - Az Y" 

M ar + Az X - Ax Z 

N ar + Ax Y - Ay X 

The total aerodynamic moment acting at the vehicle center can be expressed in 
terms of the force and moment coefficients derived from the user-supplied nonlinear 
aerodynamic modeling subroutine CCALC by defining the body axis forces in terms of 
stability axis force coefficients: 

X = qS(-C D cos a + C L sin a) 

Y = gSCy 

Z = qS(-C D sin a - Cl cos a) 


M = 


Substituting these equations into the definition of the total aerodynamic moment 
equation and applying the definitions of the total aerodynamic moments. 
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L ar — qSbC. 
M ar = q ScC m ar 


N ar ~ q s bC nar 


expressions for total aerodynamic moment coefficients corrected to the vehicle center 
of gravity can be derived as follows: 


Ay Az 

C £ = C £ + b“ ( ~ C ° sin a " C L cos <*) - ]— C Y 
ar 


Ax 

c m = c m ar + — (-C D cos a + C L sin a) - — — (-c D sin a - C L cos a) 


Az 

c 


c n = c n ar + Cy - ^ (-Cd cos a + Cl sin a) 


These calculations are normally performed within LINEAR in the subroutine 
CGCALC. However , if the user selects, the calculation can be performed within 
the user-supplied aerodynamic model, CCALC* 
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APPENDIX B: ENGINE TORQUE AND GYROSCOPIC EFFECTS MODEL 


Torque and gyroscopic effects represent (after thrust) the main contributions 
of the engines to the aircraft dynamics. The torque effects arise due to thrust 
vectors not acting at the vehicle center of gravity. The gyroscopic effects are a 
consequence of the interaction of the rotating mass of the engine and the vehicle 
dynamics. These effects can be either major or virtually negligible, depending on 
the vehicle. 

The torque effects can be modeled by considering the thrust of an engine, F p , 

where the thrust vector is aligned with the local x axis of the engine acting at 
some point Ar from the center of gravity of the vehicle, as shown in figure 4. 



z 

Figure 4. Def inition of location of engine 
center of mass ( CM relative to vehicle center 
of gravity . 


The thrust vector for the ith engine, F p ^, can be defined as 

P P yi F PzJ T 

where Fp x . / F py. / and F p z . are t ^ ie com P on ents of thrust in the x, y, and z body 

axes, respectively. From figures 5 and 6 it can be seen that the following rela- 
tionships hold: 

F p Xi = | F P1 | cos cos Si 
F P yi = I F Pi I cos e i sin ^ 

F P Zi = - I F Pi I sin e i 
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the origin of the ith thrust vector. 

The torque due to offset from the center of gravity of the ith engine, Ax 0 ^, is 
then given by 

Ax Ci = Ar^ x F p 

Thus , 


Ato. = 

i 


The total torque due to engines offset from the center of gravity of the vehicle. 
To, is given by 
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= £ At c 


i=1 


n 

= I 

i=1 


(Ari 


AF Pi> 


where n is the number of engines. 

For the case of vectored thrust, the equations for torque produced at the 
vehicle center of gravity from the ith engine, At 0 are somewhat more complicated. 

Figure 7 schematically represents an engine with thrust vectoring whose center of 
gravity is located at Ar^ relative to the vehicle's center of gravity. 



Figure 7. Detailed definition of 
thrust-vectoring parameters • 


The thrust is assummed to act at -Axrji in the local (engine) x axis, with the engine 

center of gravity being the origin of this local coordinate system. The thrust is 
also assumed to be vectored at angles hi and £i relative to the local coordinate 

axes, with hi being the angle from the thrust vector to the engine x-y plane and 
the angle from the projection of the thrust vector onto the engine x-y plane to the 
local x axis. Thus, letting F PXi , Fp^ , and Fp z> represent the x, y, and z thrust 

components in the local engine coordinate system, respectively, these terms can be 
defined in terms of the total thrust for the ith engine, Fp^, and the angles hi 
and Ci as 

= | F pi | cos Tii cos c i 
= j Fp^ | cos hi sin 
= - | F pi | sin Hi 


P Xi 

P Yi 


P Zi 
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where 


L 


To transform this equation from the ith engine 
the transformation matrix 


'cos cos 
cos Gj_ cos 


-sin 


-sin £ i 
cos 
0 


axis system to the body 


sin cos £i 
sin e-L cos £i 


axi 


O 0^0 UCill j 


cos £i 


is used. 

The resultant 


force in body axis 


is coordinates 


is 


so that 



F P Xi _ 


COS £j_ COS 5i 

-sin ?i 

sin £i cos £i 



F P1 = 

Fp Yi 

= 

cos £ sin 5i 

cos ?i 

sin £i cos £i 


F Pn 


- Fpz i. 


-sin £ 

0 

cos £i 



F Pxi = 

1 F Pi 

| (cos m cos COS £ j_ COS Ci 

- COS Tli sin ?i 

sin Ki 



- 

sin rii sin e i 

cos 5i) 




F P yi = 

1 F Pi 

| (cos tli cos ?j_ cos £ j_ sin 5i 

+ cos hi sin £i cos £i 



- 

sin rii s ^ n e i 

cos 5i) 




F P Z i ~ 

1 F Pi 

<- 

cos rii cos Ci 

sin - sin 

hi cos e i) 




The moment arm through which the vectored thrust acts is 

Ax^ - Axrp^ cos £j_ cos £i 
Ar| = | Ayj_ - Ax^ cos £i sin Ci 
Az^ + Ax«r^ sin 

and the total torque due to thrust vectoring is 

n n 

T o = X Ax °i = Z (Ar i x AF Pi> 
i=1 i=1 
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The gyroscopic effects due to the interaction of the rotating mass of the engines 
and the vehicle dynamics can be derived from the equation for the rate of change of 
angular momentum, 

T gi = ^e^ + <u x I ei w ei 

where is the gyroscopic moment produced by the ith engine, h e ^ the rate of 

change of the angular momentum of the ith engine, I Gi the inertia tensor for the ith 
engine, ou e ^ the rotational velocity of the ith engine, and w the total rotational 
velocity of the vehicle given by 


0) = [ p q r ] T 

If it is assumed that the angular momentum of the engine is constant, then 


h e . = 0 
e x 


and the equation simplifies to 


T gj_ “ w A - L e ; w e 


= uj x I 0 . au 


Two terms, I e ^ and (o e ^ , remain to be defined in the equation for the gyroscopic 

moment produced by the ith engine. Once again, simplifing assumptions are made. It 
is assumed that the inertia tensor of the engine contains a single nonzero entry. 


0 

L 0 


0 

0 

0 


0‘ 

0 

0J 


where I x is the inertia tensor of the ith engine at the location of the ith engine, 
e i 

oriented with the local x axis coincident with the rotational velocity of the engine. 
The rotational velocity of the engine has components in the x, y, and z body axes 
(p e ^, and r e^t respectively) so that 


We i ” \^ e x ^ e i re i] 


where 


P ei - 

1 w ei | 

COS €$_ COS 

5i 

q ei = 

1 “ ei | 

cos £-j_ sin 

Si 

r «i - 

- | w e 

sin ei 
i 1 A 



with u) e . being the total angular velocity of the ith engine. 


61 


The engine inertia tensor must be defined in an axis system oriented consist- 
ent with the vehicle body axis system. This is done in two steps. These steps 
involve rotating the engine inertia tensor into a coordinate system orthogonal to 
the aircraft body axis system. First, the ith engine inertia tensor is rotated 
through an angle £^ about the local y axis so that the new inertia tensor is ori- 
ented with its local x-y body axis plane parallel to the x-y body axis plane of the 
vehicle. The second step requires a rotation through an angle about the local 

z axis so that the local x, y, and z axes are orthogonal to the x, y, and z body 
axes of the vehicle. As determined by Gainer and Sherwood (1972) and Thelander 
(1965), this rotation is a similarity transformation that yields a new inertia 
tensor I e ^, such that 

Ie i = R ?i R ei I Xe i R e i R C i 

where and R e are axis transformation matrices that perform the previously 
described rotations through £ and e, respectively. These matrices are given as 


Re * = 


cos 

0 

-sin £4 


0 sin £ i 

1 0 

0 COS £j_ 


R Ci = 


cos £i 
sin 

0 


-sin 5 i 0 

cos £i 0 

0 0 


so that 



COS 

COS £ j_ 

-sin 

cos 

5i sin 

^i 

Hi R Ei = 

sin 

E, i COS £ i 

cos £i 

sin 

Ki sin 

^i 


- 

-sin £ i 

0 


COS £i 



Because 


and the matrices are unitary, 


-1 T 

R £ i - Re i 


-1 T 

Hi = H 


-1 -1 


T T 


HiHi = HjHi - (HiHi) 
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cos 5 i cos sin cos -sin £j 


( r e i R C i) T = “ sin 


cos 


Therefore, 


cos 5 i sin Si sin Ci sin cos 


cos 2 cos 2 C i cos 2 cos ?i sin ?i -cos sin e± cos £jl 


e i = Ix e . cos 2 e 


i cos 5 i sin sin 2 cos 2 -sin cos sin 


|^-cos ei sin ei cos £i 
The angular momentum of the ith 




with 


-sin e± cos e± sin £i 
th engine, h e ^, can now be expressed 

i . to a . = r*hfik . hp ■ hp ■ ”1 

i e i i e i 1 e i 2 e i 3 j 


sin 2 Ci 


e i 1 Pe i[ Ix ei 


! i[ Ix ei COs2 e i cos2 ?i ] 

q ei j^I Xe . cos 2 cos sin 5iJ 

r e .T-i x cos sin e± cos 2 

L L e i J 


cos sin e 


:i cos 2 


e i2 = Pei^Xg. cos e i cos ^i sin ?i] 

qe i[ Ix ei sin2 cos2 Ei ] 

r e . r-i x sin cos 5i sin 

i e-; 


,[ 1 x ei s 111 *" s i cos 2 c-i 

[ -I x sin cos Si 

e i 


i sin 5iJ 


le i 3 “ Pe i " Ix ei COS £i 


■ j_ sin £i cos 2 


+ q ei [-Ix ei sin e i cos sin Si] 

+ r «i[ Ix ei Sin2 Ei ] 


sin 


gyroscopic moment induced by the ith engine, Tg , can be expanded 


r q he i3 - rh e 


Tg. = rh eii - ph e . 3 


Ph eil “ qh ei2 
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and the total moment induced by gyroscopic interaction of the vehicle dynamics 
and the rotating engine components is 


n 

T g = X T< ?i 

i = 1 

Engine torque and gyroscopic effects are modeled within the subroutine ENGINE 
using information provided by the user from the engine modeling subroutine IFENGN. 
These effects are calculated as incremental moments and are included directly in 
the equations of motion for both analysis point definition and derivation of the 
linearized system matrices. 
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APPENDIX C: STATE VARIABLE NAMES RECOGNIZED BY LINEAR 


This appendix lists the alphanumeric descriptors specifying state variables that 
are recognized by LINEAR. In the input file, the field containing these descriptors 
uses a 54A format, and all characters are left-justified. The input alphanumeric 
descriptor specified by the user serves both to identify the state variable selected 
by the user within the program itself and to identify state variables on the printed 
output of LINEAR, as described in the Output Files section. 


State variable 

Units 

Symbol 

Alphanumeric 

descriptor 

Roll rate 

rad/sec 

P 

p 

ROLL RATE 

Pitch rate 

rad/sec 

q 

Q 

PITCH RATE 

Yaw rate 

rad/sec 

r 

R 

YAW RATE 

Velocity 

ft/sec 

V 

V 

VELOCITY 

VEL 

VTOT 

Angle of 
attack 

rad 

a 

ALP 

ALPHA 

ANGLE OF ATTACK 

Sideslip angle 

rad 

3 

BTA 

BETA 

SIDESLIP ANGLE 
ANGLE OF SIDESLIP 

Pitch angle 

rad 

0 

THA 

THETA 

PITCH ATTITUDE 

Heading 

rad 

♦ 

PS I 

HEADING 
HEADING ANGLE 

Roll angle 

rad 


PHI 

ROLL ATTITUDE 
BANK ANGLE 

Altitude 

length 

h 

H 

ALTITUDE 

Displacement 

north 

length 

X 

X 

Displacement 

east 

length 

y 

Y 
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APPENDIX D: OBSERVATION VARIABLE NAMES RECOGNIZED BY LINEAR 


This appendix lists all observation variable names recognized by LINEAR except 
for state and control variable names. If state variables are specified as elements 
in the observation vector, the alphanumeric descriptor must correspond to the names 
defined in appendix C. When control variables are to be included in the observation 
vector these variables must be identified exactly as they were specified by the user. 

The input file is formatted 5A4 with the alphanumeric data left- justified . The 
floating-point fields are used to define sensor locations not at the center of grav- 
ity. The input name specified by the user for an observation variable serves both 
to identify the observation variable selected within the program itself and to iden- 
tify observation variables on the printed output of LINEAR. 

An asterisk preceding the variable name indicates measurements at some point 
other than the vehicle center of gravity. The program LINEAR uses the quantities 
defined in the first three floating-point fields as definitions of the location 
of the sensor with respect to the vehicle center of gravity. The three parameters 
define the x body, y body, and z body location, in that order, of the sensor. 

These offsets from the vehicle center of gravity are defined in units of length. 


Observation variable 

Units 

Symbol 

Alphanumeric descriptor 

Derivatives of state variables 

Roll acceleration 

rad/sec 2 

• 

P 

PDOT 

ROLL ACCELERATION 

Pitch acceleration 

rad/sec 2 

• 

q 

QDOT 

PITCH ACCELERATION 

Yaw acceleration 

rad/sec 2 

• 

r 

RDOT 

YAW ACCELERATION 

Velocity rate 

length/sec 2 

• 

V 

VDOT 

VELOCITY RATE 

Angle-of-attack rate 

rad/sec 

• 

a 

ALPDOT 
ALPHA DOT 
ALPHADOT 

Angle-of-sideslip rate 

rad/sec 

• 

e 

BTADOT 
BETA DOT 
BETADOT 

Pitch attitude rate 

rad/sec 

• 

0 

THADOT 
THETA DOT 

Heading rate 

rad/sec 

4 1 

PSIDOT 
PS I DOT 
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Observation variable 

Units 

Symbol 

Alphanumeric 

descriptor 

Derivatives of state 

variables 

(continued ) 


Roll attitude rate 

rad/sec 

• 

<t> 

PHIDOT 
PHI DOT 


Altitude rate 

length/sec 

• 

h 

HDOT 

ALTITUDE RATE 

Velocity north 

length/sec 

• 

X 

XDOT 


Velocity east 

length/sec 

• 

y 

YDOT 



Accelerations 


x body axis 
acceleration 

g 

a x 

AX 

LONGITUDINAL ACCEL 
X-AXIS ACCELERATION 
X AXIS ACCELERATION 
X-BODY AXIS ACCEL 
X BODY AXIS ACCEL 

y body axis 
acceleration 

g 

a y 

AY 

Y-AXIS ACCELERATION 

Y AXIS ACCELERATION 
Y-BODY AXIS ACCEL 

Y BODY AXIS ACCEL 
LATERAL ACCELERATION 
LAT ACCEL 

LATERAL ACCEL 

z body axis 
acceleration 

g 

a z 

AZ 

Z-BODY AXIS ACCEL 
Z BODY AXIS ACCEL 

x body axis acceler- 
ometer at vehicle 
center of gravity 

g 

a nx 

ANX 

X-AXIS ACCELEROMETER 
X AXIS ACCELEROMETER 

y body axis acceler- 
ometer at vehicle 
center of gravity 

g 

a ny 

ANY 

Y-AXIS ACCELEROMETER 
Y AXIS ACCELEROMETER 


Observation variable 

Units 

Symbol 

Alphanumeric descriptor 


Accelerations 

(continued ) 


z body axis acceler- 
ometer at vehicle 
center of gravity 

g 

a nz 

ANZ 

Z-AXIS ACCELEROMETER 
Z AXIS ACCELEROMETER 

Normal acceleration 

g 

a„ 

n 

AN 

NORMAL ACCELERATION 
NORMAL ACCEL 
GS 
G'S 

*x body axis accel- 
erometer not at 
vehicle center of 
gravity 

g 

a nx, i 

AX, I 
ANX, I 

*y body axis accel- 
erometer not at 
vehicle center of 
gravity 

g 

a ny,i 

AY , I 
ANY, I 

*z body axis accel- 
erometer not at 
vehicle center of 
gravity 

g 

a nz ,i 

AZ , I 
ANZ , I 

♦Normal accelerometer 
not at vehicle 
center of gravity 

g 

a n,i 

AN, I 

Load factor 

(Dimension- 

less) 

n 

N 

LOAD FACTOR 

Air data parameters 

Speed of sound 

length/sec 

a 

A 

SPEED OF SOUND 

a Reynolds number 

(Dimension- 

less 

Re 

RE 

REYNOLDS NUMBER 

Reynolds number 
per unit length 

length *" 1 

Re ' 

RE PRIME 
R/LENGTH 
R/FEET 

R/UNIT LENGTH 
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a Reynolds number is defined in terms of an arbitrary unit of length that is 
input by the user. This length is input using the first floating-point field 

however, if no value is input, c is used as the default value. 



Observation variable 


Units 


Symbol Alphanumeric descriptor 



Air data parameters 

(continued ) 

Mach number 

(Dimension- 

M 

M 


less ) 


MACH 

AMCH 

Dynamic pressure 

force/ 

q 

QBAR 


length 2 


DYNAMIC PRESSURE 

Impact pressure 

force/ 

<?c 

QC 


length 2 


IMPACT PRESSURE 
DIFFERENTIAL PRESSURE 

Ambient pressure 

force/ 

Pa 

PA 


length 2 


STATIC PRESSURE 
FREES TREAM PRESSURE 

Impact/ambient 

(Dimension- 

^c/Pa 

QC/PA 

pressure ratio 

less ) 


QC/ P 

Total pressure 

force/ 

Pt 

PT 


length 2 


TOTAL PRESSURE 

Temperature 

K 

T 

TEMP 

TEMPERATURE 
FREESTREAM TEMPERATURE 

Total temperature 

K 

T t 

TOTAL TEMPERATURE 

Equivalent airspeed 

knots 

V e 

VEAS 

EQUIVALENT AIRSPEED 
KEAS 

Calibrated airspeed 

knots 

v c 

VCAS 

CALIBRATED AIRSPEED 
KCAS 


Flightpath-related 

parameters 

Flightpath angle 

rad 

Y 

GAM 

GAMMA 

FLIGHT PATH ANGLE 
GLIDE PATH ANGLE 
GLIDE SLOPE 

Flightpath 

g 

fpa 

FPA 

acceleration 



FLIGHT PATH ACCEL 


Observation variable 

Units 

Symbol 

Alphanumeric descriptor 

Flightpath-related 

parameters ( continued ) 

Flightpath angle 
rate 

rad/sec 

• 

Y 

GAMMA DOT 
GAMMADOT 

Vertical acceler- 
ation 

length/sec 

2 h 

VERTICAL ACCELERATION 
HDOTDOT 
H -DOT -DOT 
HDOT— DOT 

Scaled altitude 
rate 

length/sec 

h/57.3 

H -DOT/5 7. 3 
HDOT / 57.3 


Energy-: 

related terms 


Specific energy 

length 

E S 

ES 

E-SUB-S 

SPECIFIC ENERGY 

Specific power 

length/sec 


PS 

P-SUB-S 

SPECIFIC POWER 
SPECIFIC THRUST 


Force 

parameters 


Lift force 

force 

L 

LIFT 

Drag force 

force 

D 

DRAG 

Normal force 

force 

N 

NORMAL FORCE 

Axial force 

force 

A 

AXIAL FORCE 

Body axis parameters 

x body axis velocity 

length/sec 

u 

UB 

X-BODY AXIS VELOCITY 
X BODY AXIS VELOCITY 
X-BODY AXIS VEL 
X BODY AXIS VEL 
U-BODY 
U BODY 



Observation variable 

Units 

Symbol 

Alphanumeric descriptor 

Body axis parameters (continued) 

y body axis velocity 

length/sec 

V 

VB 

Y-BODY AXIS VELOCITY 

Y BODY AXIS VELOCITY 
Y-BODY AXIS VEL 

Y BODY AXIS VEL 
V-BODY 

Y BODY 

z body axis velocity 

length/sec 

w 

WB 

Z-BODY AXIS VELOCITY 
Z BODY AXIS VELOCITY 
Z-BODY AXIS VEL 
Z BODY AXIS VEL 
W-BODY 
W-BODY 

Rate of change of 
velocity in x body 
axis 

length/sec 2 

• 

u 

UBDOT 
UB DOT 

Rate of change of 
velocity in y body 
axis 

length/sec 2 

• 

V 

VBDOT 
VB DOT 

Rate of change of 
velocity in z body 
axis 

length/sec 2 

• 

w 

WBDOT 
WB DOT 

Miscellaneous measurements not 

at vehicle 

center of gravity 

*Angle of attack not 
at vehicle center 
of gravity 

rad 

a ,i 

ALPHA, I 

ALPHA INSTRUMENT 
AOA INSTRUMENT 

*Angle of sideslip 
not at vehicle 
center of gravity 

rad 

®.i 

BETA, I 

BETA INSTRUMENT 
SIDESLIP INSTRUMENT 

* Altitude instrument 
not at vehicle 
center of gravity 

length 

h ,i 

H,I 

ALTITUDE INSTRUMENT 

^Altitude rate instru- 
ment not at vehicle 
center of gravity 

length/sec 


HDOT , I 

• 


Observation variable 

Units 

Symbol 

Alphanumeric descriptor 

Other miscellaneous parameters 

Vehicle total angular 
momentum 

mass-length 2 / 

sec 2 

T 

ANGULAR MOMENTUM 
ANG MOMENTUM 

Stability axis roll 
rate 

rad/sec 

Ps 

STAB AXIS ROLL RATE 

Stability axis pitch 
rate 

rad/sec 

<Is 

STAB AXIS PITCH RATE 

Stability axis yaw 
rate 

rad/sec 

r s 

STAB AXIS YAW RATE 




♦ 
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APPENDIX E: ANALYSIS POINT DEFINITION IDENTIFIERS 


Analysis point definition options are selected using alphanumeric descriptors. 
The first record read for each analysis case contains these descriptors. All these 
descriptors are read using a 5A4 format. The following table associates the analy- 
sis point definition options with their alphanumeric descriptors. 


Analysis point definition option 

Alphanumeric descriptor 

Un trimmed 

UNTRIMMED 
NO TRIM 
NONE 
NOTRIM 

Straight-and-level 

STRAIGHT AND LEVEL 
WINGS LEVEL 
LEVEL FLIGHT 

Pushover-pullup 

PUSHOVER/PULLUP 

PULLUP 

PUSHOVER 

PUSHOVER AND PULLUP 
PUSHOVER / PULLUP 
PUSHOVER / PULL-UP 
PUSHOVER PULLUP 
PUSH OVER 
PULL UP 

Level turn 

LEVEL TURN 
WINDUP TURN 

Thrust-stabilized turn 

THRUST STABILIZED TURN 
THRUST LIMITED TURN 
FIXED THROTTLE TURN 
FIXED THRUST TURN 

Beta 

BETA 


SIDESLIP 

Specific power 

SPECIFIC POWER 
PS 

P-SUB-S 


Each of these analysis point definitions except the un trimmed, beta, and spe- 
cific power options has two suboptions associated with it. The suboptions are 
requested using alphanumeric descriptors read using an A4 format. These suboptions 
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are defined in the Analysis Point Definition section. The following table defines 
these suboptions and the alphanumeric descriptors associated with each. 


Analysis point definition suboption 

Alphanumeric descriptor 

Straight-and-level 


Alpha-trim 

ALP 

ALPH 

ALPHA 

Mach- trim 
Pushover-pullup 

MACH 

AMCH 

Alpha-trim 

ALP 

ALPH 

ALPHA 

Load-factor-trim 

LOAD 
GS 
G 1 S 
AN 

Level turn 


Alpha-trim 

ALP 

ALPH 

ALPHA 

Load-f actor- trim 

LOAD 
GS 
G *S 
AN 

Thrust-stabilized turn 


Alpha- trim 

ALP 

ALPH 

ALPHA 

Load-f actor- trim 

LOAD 
GS 
G * S 
AN 


I 


74 


APPENDIX F: EXAMPLE INPUT FILE 


The following listing is an example of an input file to LINEAR. This file was 
used with the example subroutines listed in appendix I to generate the analysis and 
printer output files listed in appendixes G and H, respectively. 


LINEARIZER TEST AND DEMONSTRATION CASES 
USER’S GUIDE 

6. 080000E+02 4.280000E+01 1.595000E+-1 4.500000E+04 
2. 870000E+04 1.651000E+-5 1.879000E+05-5. 200000E+02 0. 0. 

0. 0. 0. CCALC WILL CALCULATE CG CORRECTIONS 

-1.000000E+01 4. OOOOOOE+Ol 

4STAN RADI 

ALPHA 

Q 

THETA 

VEL 

3 

ELEVATOR 5 

THROTTLE 12 

SPEED BRAKE 10 

2STAN 
AN 
AY 

-2. 900E+00 5, 430E+00-4.000E+00 4.000E+00-3.250E+00 3. 250E+00-1 .OOOE+OO 1.000E+00 
0 ADDITIONAL SURFACES 


WINDUP TURN 

ALPHA 

H 

20000.0 

MACH 

0.90 

AN 

3.00 

BETA 

0.0 

NEXT 

LEVEL FLIGHT 

ALPHA 

H 

20000.0 

MACH 

0.9 

GAMMA 

10.0 

END 



This input file is for a case called (record 1) LINEARIZER TEST AND DEMONSTRA- 
TION CASES r and all input data are on logical device unit 1, signified by the second 
record being blank. The project title is USER'S GUIDE (record 3). Record 4 
specifies the mass and geometric properties of the vehicle as 


S = 608 ft 2 
b = 42.8 ft 


c = 15.95 ft 
w = 45,000 lb 

Record 5 defines the moments and products of intertia of the vehicle (all in 
units of slug-ft^) as 


I 
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I x = 28,700 
I y = 165,100 
I z = 187,900 
^xz = 520 
I X y = 0 
lyz = 0 

Record 6 defines the location of the vehicle center of gravity to be coincident with 
the aerodynamic reference point of the nonlinear aerodynamic model by setting 

Ax = 0 

Ay = 0 

Az = 0 

Record 6 also specifies that LINEAR should not use its internal model to make cor- 
rections for the offset in the vehicle center of gravity from the aerodynamic ref- 
erence point, either because the aerodynamic model includes such corrections or 
because none are to be made. 

Record 7 defines the angle-of-attack range of the aerodynamic model. 

Record 8 specifies that there will be four state variables in the output, 
that the output formulation of the state equation will be in the standard form 

(x = Ax + Bu), and that the output for the nondimensional stability and control 
derivatives with respect to angle of attack and angle of sideslip should be scaled 
as reciprocal radians. 

The next four records (9 to 12) define the output formulation of the state 
vector to be 


a 


x = 


q 

e 


V 


Record 13 specifies that the output model will have three parameters in the con- 
trol vector. The following three records (14 to 16) specify that 


u 


elevator 
throttle 
speed brake 


and that elevator, throttle, and speed brake are located in DC(5), DC(12), and DC(10) 
of the /CONTROL/ common block. 
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Record 17 specifies that two observation variables will be output and that the 
observation equation will be in the standard form, 

y = Hx + Fu 

The next two records (18 and 19) define the elements of the output vector to be 



Record 20 specifies the ranges for the trim parameters, DES, DAS, DRS , and 
THRSTX, used to trim the longitudinal, lateral, and directional axes and thrust, 
respectively. The ranges for these parameters are defined by record 20 to be 

-2.9 < DES < 5.43 
-4.0 < DAS < 4.0 
-3.25 < DRS < 3.25 
-1 .0 < THRSTX < 1.0 

The first three parameters essentially represent stick and rudder limits and are 
so specified because of the implementation of the subroutine UCNTRL (discussed in 
app. I). The thrust trim parameter is used, again because of the implementation 
of UCNTRL, to schedule speed brake when THRSTX < 0 and to command thrust when 
THRSTX > 0. 

Record 21 specifies that no additional control surfaces are to be set. 

The next seven records (22 to 28) define an analysis point option. These records 
request a level turn trim option at 


h = 20,000 ft 
M = 0.9 
a n = 3.0 g 
5 = 0 

The second record of this set (record 23) indicates which level turn suboption is 
requested. The alphanumeric descriptor ALPHA indicates that angle of attack is to 
be varied until the specified 3.0-g turn is achieved. The final record of this anal- 
ysis point option definition set contains the key word NEXT to indicate both an 
end to the current analysis point option definition and that another analysis point 
option definition follows. 

The final six records (records 29 to 34) define a straight-and-level analysis point 
option at 


h = 20,000 ft 
M = 0.9 
Y = 10.0° 
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The second record of this set (record 30) identifies the M Alpha-trim" suboption in 
which angle of attack is varied until trim is achieved at the specified condition. 
The final record of this set contains the key word END to indicate the termination 
of the current analysis point definition as well as the termination of input cases. 


78 



APPENDIX G: EXAMPLE OUTPUT ANALYSIS FILE 


The following listing is an example analysis file output on unit 15. This file 
was produced using the example input file listed in appendix F and the example user- 
supplied subroutines listed in appendix I. 


LINEARIZER TEST AND DEMONSTRATION CASES 
USER'S GUIDE 

TEST CASE # ************************************************* i 

X - DIMENSION =4 U - DIMENSION =3 Y - DIMENSION = 2 
STATE EQUATION FORMULATION: STANDARD 

OBSERVATION EQUATION FORMULATION: STANDARD 


STATE VARIABLES 


ALPHA 

= 

0.465695D-01 

RADIANS 

Q 

= 

0. 921683D-01 

RADIANS/SECOND 

THETA 

= 

0. 159885D -01 

RADIANS 

VEL 

= 

0. 933232D+03 

FEET/SECOND 

CONTROL VARIABLES 

ELEVATOR 


0. 538044D -01 


THROTTLE 

= 

0.214105D+00 


SPEED BRAKE 

= 

O.OOOOOOD+OO 



DYNAMIC INTERACTION VARIABLES 


X-BODY AXIS FORCE 

= O.OOOOOOD+OO 

POUNDS 

Y-BODY AXIS FORCE 

= O.OOOOOOD+OO 

POUNDS 

Z-BODY AXIS FORCE 

= O.OOOOOOD+OO 

POUNDS 

PITCHING MOMENT 

= O.OOOOOOD+OO 

FOOT-POUNDS 

ROLLING MOMENT 

= O.OOOOOOD+OO 

FOOT-POUNDS 

YA'JING MOMENT 

= O.OOOOOOD+OO 

FOOT-POUNDS 


OBSERVATION VARIABLES 

0. 300163D+01 GS 
0.941435D+00 GS 


AN 

AY 


A-MATRIX FOR: DX/DT = A * X + B*U+D*V 

-0. 121436D+01 0. 100000D+01 0. 1367 56D -02-0. 121 605D -03 
-0. 147423D+01 -0.221451D+01-0. 450462D-02 0.294019D-03 
O.OOOOOOD+OO 0. 331812D+00 O.OOOOOOD+OO O.OOOOOOD+OO 
-0. 790853D+02 0.000000D+O0-0.320822D+02-0. 157297D-01 


B-MATRIX FOR: DX/DT = A * X + B*U+D # V 

-0. 141961D+00-0. 1649480 -02-0. 92 8933D-02 
-0. 220778D+02 0. 543324D -02-0. 135074D+02 
O.OOOOOOD+OO O.OOOOOOD+OO O.OOOOOOD+OO 
-0. 105186D+02 0. 34281 7D+02-0. 155832D+02 


D-MATRIX FOR: DX/DT=A*X+B*U+D*V 


-0. 343642D -07 
0. 113192D-06 
O.OOOOOOD+OO 
0. 714203D -03 


O.OOOOOOD+OO 0. 737378D-06 O.OOOOOOD+OO O.OOOOOOD+OO 
0 . 000000D +00 -0 . 242 885D -05 0.605694D-05 O.OOOOOOD+OO 
O.OOOOOOD+OO O.OOOOOOD+OO O.OOOOOOD+OO O.OOOOOOD+OO 
0. 3984 92D -06 0.332842D-04 O.OOOOOOD+OO O.OOOOOOD+OO 
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0. OOOOOOD+OO 
0. 000000D+00 
0. OOOOOOD+OO 
0. OOOOOOD+OO 


H-MATRIX FOR: Y = H * X + F * U + E • V 

0.351752D+02-0. 191245D-05 0.150046D-02 0.640771D-02 
0. OOOOOOD+OO 0* OOOOOOD+OO -0, 150534D -01 0.274248D-09 


F -MATRIX FOR: Y = H • X + F # U+E*V 

0. 41284 5D+01 -0. 180978D -02 0. 291699D+00 
0. OOOOOOD+OO 0. OOOOOOD+OO 0. OOOOOOD+OO 


E-MATRIX FOR: Y = H*X+F*U+E*V 

-0. 377037D -07 0. OOOOOOD+OO -0, 2141 32D -04 0. OOOOOOD+OO 0. OOOOOOD+OO 
0. OOOOOOD+OO 0, 222222D -04 0. OOOOOOD+OO 0. OOOOOOD+OO 0. OOOOOOD+OO 
0* OOOOOOD+OO 
0. OOOOOOD+OO 


TEST CASE ************************************************ 2 

X - DIMENSION =4 U - DIMENSION =3 Y - DIMENSION = 2 
STATE EQUATION FORMULATION: STANDARD 

OBSERVATION EQUATION FORMULATION: STANDARD 


STATE VARIABLES 

ALPHA 

= 

-. 126650D -01 

RADIANS 

Q 

= 

0. OOOOOOD+OO 

RADIANS/SECOND 

THETA 


0. 16186 8D +00 

RADIANS 

VEL 

= 

0. 933232D+03 

FEET/ SECOND 

CONTROL VARIABLES 

ELEVATOR 

= 

0. 637734D -01 


THROTTLE 

= 

0. 225092D+00 


SPEED BRAKE 

= 

0. OOOOOOD+OO 



DYNAMIC INTERACTION VARIABLES 


X-BODY AXIS FORCE 

= 0. OOOOOOD+OO 

POUNDS 

Y-BODY AXIS FORCE 

= 0. OOOOOOD+OO 

POUNDS 

Z-BODY AXIS FORCE 

= 0. OOOOOOD+OO 

POUNDS 

PITCHING MOMENT 

= 0. OOOOOOD+OO 

FOOT-POUNDS 

ROLLING MOMENT 

= 0. OOOOOOD+OO 

FOOT-POUNDS 

YAWING MOMENT 

= 0. OOOOOOD+OO 

FOOT-POUNDS 


OBSERVATION 

VARIABLES 


AN 

= 0. 985228D+O0 

GS 

AY 

= 0. OOOOOOD+OO 

GS 

A-MATRIX FOR 

: DX / DT = A * X 

+ B * U + 1 


-0. 12 0 900D+01 0. 100000D+01 -0. 575730D -02-0. 701 975D -04 
-0. 149189D+01 -0. 221451D+01 0.189640D-01 0.231368D-03 
0. OOOOOOD+OO 0. 100000D+01 0. OOOOOOD+OO 0. OOOOOOD+OO 
-0. 576868D+02 0. OOOOOOD+OO -0. 31 62 BID +02 -0.4604 3 5D -02 


<b <b 


B -MATRIX FOR: DX / DT 


A * X + B * U + D * V 


. 141 961D+00 0. 448742D -03 -0. 928932D -02 
. 220778D+02 -0. 147812D -02 -0. 135074D+02 
0. OOOOOOD+OO 0. OOOOOOD +00 O.OOOOOOD+OO 
-0, 105186D+02 0. 343162D+02-0. 155832D+02 


D-MATRIX FOR: DX/DT+A*X+B*U+D*V 

0. 934880D-08 O.OOOOOOD+OO 0. 738119D-06 O.OOOOOOD+OO O.OOOOOOD+OO 
-0. 307941D -07 0. OOOOOOD+OO -0. 2431 2 9D -05 0.605694D-05 O.OOOOOOD+OO 
O.OOOOOOD+OO O.OOOOOOD+OO 0. 000000D-K)0 O.OOOOOOD+OO O.OOOOOOD+OO 
0. 71492 OD -0 3 0 . 495 82 90 -13-0 . 90 5 4 97D -05 O.OOOOOOD+OO O.OOOOOOD+OO 
O.OOOOOOD+OO 
O.OOOOOOD+OO 
O.OOOOOOD+OO 
O.OOOOOOD+OO 


H-MATRIX FOR: Y = H * X + F * U 

0. 350424D+02 -0 . 12 8333D -06 -0 . 632314D -02 0. 203434D -02 
O.OOOOOOD+OO O.OOOOOOD+OO O.OOOOOOD+OO O.OQOOOOD+OO 

F -MATRIX FOR: Y = H * X + F * U 

0.411323D+0 1 0.492 84 5D -03 0.263288D+00 
O.OOOOOOD+OO O.OOOOOOD+OO O.OOOOOOD+OO 


E -MATRIX FOR: Y + H*X+F*U+E*V 

0. 102676D-07 0. OOOOOOD+OO -0. 21 411 6D -04 O.OOOOOOD+OO O.OOOOOOD+OO 
O.OOOOOOD+OO 0.222222D -04 O.OOOOOOD+OO O.OOOOOOD+OO O.OOOOOOD+OO 
O.OOOOOOD+OO 
O.OOOOOOD+OO 


APPENDIX H: EXAMPLE PRINTER OUTPUT FILES 


The following listings are the printer output files generated by LINEAR using 
the example input file in appendix F and the example user supplied subroutines 
listed in appendix I. 


Example printer output file 1 (unit 3) 


GEOMETRY AND MASS DATA FOR: 

LINEARIZER TEST AND DEMONSTRATION CASES 


FOR THE PROJECT: 
USER'S GUIDE 


WING 

AREA 

608.000 

(FT**2 ) 

WING 

SPAN 

42.800 

(FT) 

MEAN 

CHORD 

15.950 

(FT) 

VEHICLE WEIGHT 

45000.000 

(LB) 

IX 


28700.000 

(SLUG-FT**2) 

IY 


165100.000 

(SLUG -FT**2 ) 

IZ 


187900.000 

(SLUG-FT**2 ) 

IXZ 


-520.000 

(SLUG-FT**2 ) 

IXY 


0.000 

(SLUG-FT**2 ) 

IYZ 


0.000 

(SLUG-FT**2 ) 


VECTOR DEFINING REFERENCE POINT OF AERODYNAMIC MODEL 
WITH RESPECT TO VEHICLE CENTER OF GRAVITY: 


DELTA X 0.000 (FT) 

DELTA Y 0.000 (FT) 

DELTA Z 0.000 (FT) 

FORCE AND MOMENT COEFFICIENT CORRECTIONS DUE 
TO THE OFFSET OF THE REFERENCE POINT TO THE 
AERODYNAMIC MODEL FROM THE VEHICLE CG ARE 
CALCULATED IN CCALC. 

MINIMUM ANGLE OF ATTACK -10.000 (DEG) 

MAXIMUM ANGLE OF ATTACK 40.000 (DEG) 


PARAMETERS USED IN THE STATE VECTOR FOR: 
LINEARIZER TEST AND DEMONSTRATION CASES 

FOR THE PROJECT: 

USER'S GUIDE 

ALPHA 

Q 

THETA 

VEL 


THE STANDARD FORMULATION OF THE STATE EQUATION 
HAS BEEN SELECTED. THE FORM OF THE EQUATION IS: 

DX/DT = A*X+B*U+D*V 

SURFACES TO BE USED FOR THE CONTROL VECTOR FOR: 

LINEARIZER TEST AND DEMONSTRATION CASES 

FOR THE PROJECT: 

USER'S GUIDE 

LOCATION IN /CONTROL/ 

ELEVATOR 5 

THROTTLE 12 

SPEED BRAKE 10 
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PARAMETERS USED IN THE OBSERVATION VECTOR FOR: 
LINEAR IZER TEST AND DEMONSTRATION CASES 

FOR THE PROJECT: 

USER’S GUIDE 

AN 

AY 


THE STANDARD FORMULATION OF THE OBSERVATION EQUATION 
HAS BEEN SELECTED. THE FORM OF THE EQUATION IS: 

Y=H*X+F*U+E*V 


LIMITS FOR TRIM OUTPUT PARAMETERS: 

MINIMUM MAXIMUM 
PITCH AXIS PARAMETER -2.900 “5.430 
ROLL AXIS PARAMETER -4.000 4.000 
YAW AXIS PARAMETER -3.250 3.250 
THRUST PARAMETER -1.000 1.000 


NO ADDITIONAL SURFACES TO BE SET WERE DEFINED 


TRIM CONDITIONS FOR CASE # 1 

LINEARIZER TEST AND DEMONSTRATION CASES 

FOR THE PROJECT: 

USER’S GUIDE 

LEVEL TURN WHILE VARYING ALPHA 
TRIM ACHIEVED 


COEFFICIENT OF LIFT 


= 

0.40144 

COEFFICIENT OF DRAG 


= 

0.03058 

LIFT 

(LBS) 


134741.68807 

DRAG 

(LBS) 

= 

10265.70661 

ALTITUDE 

(FT) 

= 

20000. 

MACH 


= 

0.90000 

VELOCITY 

(FT/ SEC) 


933.23196 

EQUIVALENT AIRSPEED 

(KTS) 

= 

403.42303 

SPEED OF SOUND 

(FT/SEC) 

= 

1036.92440 

GRAVITATIONAL ACCEL 

(FT/SEC**2 ) 


32.11294 

NORMAL ACCELERATION 

(G-S) 

= 

3.00163 

LOAD FACTOR 


= 

2.99995 

DYNAMIC PRESSURE 

(LBS/FT**2) 

s 

552.05302 

DENSITY 

(SLUG/FT**3) 

= 

0.00126774 

WEIGHT (OALTITUDE ) 

(LBS) 

= 

44914.60434 

BETA 

(DEG) 

= 

0.03193 

ALPHA 

(DEG) 

= 

2.66824 

PHI 

(DEG) 

= 

70.62122 

THETA 

(DEG) 

ss 

0.91607 

ALTITUDE RATE 

(FT/ SEC) 

s 

0.00000 

GAMMA 

(DEG) 

= 

0.00000 

ROLL RATE 

(DEG/SEC) 

= 

-0.08951 

PITCH RATE 

(DEG/SEC) 

= 

5.28086 

YAW RATE 

(DEG/SEC) 

= 

1.85749 

THRUST 

(LBS) 

s 

10277.03515 

SUM OF THE SQUARES 

= 

0.00000 

TRIM PARAMETERS 

TRIM PITCH AXIS PARAMETER 


-0.66958 

TRIM ROLL AXIS PARAMETER 


-0.01526 

TRIM YAW AXIS PARAMETER 

= 

-0.02125 

TRIM THRUST PARAMETER 

= 

0.21410 



CONTROL VARIABLES 


ELEVATOR 


0.05380 


THROTTLE 

= 

0.21410 


SPEED BRAKE 

- 

0.00000 


OBSERVATION VARIABLES 

AN 

- 

3.00163339 

GS 

AY 

= 

0.94136286 

GS 


NON-DIMENS IONAL STABILITY AND CONTROL DERIVATIVES FOR CASE # 1 

LINEARIZER TEST AND DEMONSTRATION CASES 


FOR THE PROJECT: 
USER'S GUIDE 




ROLLING 

PITCHING 

YAWING 



SIDE 



MOMENT 

MOMENT 

MOMENT 

DRAG 

LIFT 

FORCE 

ZERO COEFFICIENTS 

-4. 02966D -05 

4. 22040D -02 

2. 25747D -04 

1.42882D -04 

1. 57360D -01 

5.42 882D-04 

ROLL RATE 

(RAD/SEC) 

-2.00000D -01 

0. OOOOOD+OO 

-3.37210D -02 

0. OOOOOD+OO 

0. OOOOOD+OO 

0. OOOOOD+OO 

PITCH RATE 

(RAD/SEC) 

0. OOOOOD+OO 

3. 89530D+00 

0. OOOOOD+OO 

0. OOOOOD+OO 

-1 . 72320D+01 

0. OOOOOD+OO 

YAW RATE 

(RAD/SEC) 

1. 50990D -01 

0. OOOOOD+OO 

-4.04710D -01 

0. OOOOOD+OO 

0. OOOOOD+OO 

0. OOOOOD+OO 

VELOCITY 

(FT/ SEC) 

-1.27955D-07 

-3.28739D-06 

3. 21096D -07 

0. OOOOOD+OO 

1 . 45433D -05 

0. OOOOOD+OO 

MACH NUMBER 


-1 . 32680D -04 

-3. 40878D -03 

3. 32952D -04 

0. OOOOOD+OO 

1 . 50803D -02 

0. OOOOOD+OO 

ALPHA 

(RAD) 

0. OOOOOD+OO 

-1 . 6882 OD -01 

0. OOOOOD+OO 

3. 72570D -01 

4. 87060D+00 

0. OOOOOD+OO 

BETA 

(RAD) 

-1.33450D-01 

0. OOOOOD+OO 

1 . 29960D -01 

0. OOOOOD+OO 

0. OOOOOD+OO 

-9. 74030D -01 

ALTITUDE 

(FT) 

0. OOOOOD+OO 

0. OOOOOD+OO 

0. OOOOOD+OO 

0. OOOOOD+OO 

0. OOOOOD+OO 

0. OOOOOD+OO 

ALPHA DOT 

(RAD/SEC) 

0. OOOOOD+OO 

-1.18870D+O1 

0. OOOOOD+OO 

0. OOOOOD+OO 

1. 72320D+01 

0. OOOOOD+OO 

BETA DOT 

(RAD/SEC) 

0. OOOOOD+OO 

0. OOOOOD+OO 

0. OOOOOD+OO 

0. OOOOOD+OO 

0. OOOOOD+OO 

0. OOOOOD+OO 

ELEVATOR 


0. OOOOOD+OO 

-6* 95280D -01 

0. OOOOOD+OO 

4. 38310D-02 

5. 72960D -01 

0. OOOOOD+OO 

THROTTLE 


0. OOOOOD+OO 

0. OOOOOD+OO 

0. OOOOOD+OO 

0. OOOOOD+OO 

0. OOOOOD+OO 

0. OOOOOD+OO 

SPEED BRAKE 


0. OOOOOD+OO 

-4. 17500D -01 

0. OOOOOD+OO 

6. 49350D -02 

3. 74 92 OD -02 

0. OOOOOD+OO 


VEHICLE STATIC MARGIN IS 3*5% MEAN AERODYNAMIC CHORD STABLE 

AT THIS FLIGHT CONDITION. 


MATRIX A USING THE FORMULATION OF THE STATE EQUATION: 
DX/DT = A * X + B * U + D * V 
FOR CASE # 1 

LINEARIZER TEST AND DEMONSTRATION CASES 

FOR THE PROJECT: 

USER'S GUIDE 


DERIVATIVES 

WRT 

TIME OF: 

ALPHA 

Q 

THETA 

VEL 

ALPHA 

-0. 121436D+01 

0. 100000D+01 

0. 136756D-02 

-0. 121605D-03 

Q 

-0. 147423D+01 

-0. 221451D+01 

-0.450462D-02 

0. 29401 90-03 

THETA 

O.OOOOOOD+OO 

0. 331812D+O0 

O.OOOOOOD+OO 

O.OOOOOOD+OO 

VEL 

-0. 790853D+O2 

O.OOOOOOD+OO 

-0. 320822D+02 

-0. 157297D -01 


MATRIX B USING THE FORMULATION OF THE STATE EQUATION: 
DX/DT = A * X + B * U + D * V 
FOR CASE # 1 

LINEARIZER TEST AND DEMONSTRATION CASES 

FOR THE PROJECT: 

USER'S GUIDE 
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ELEVATOR THROTTLE SPEED BRAKE 

DERIVATIVES 

WRT 

TIME OF: 

ALPHA -0. 141961D+00 -0. 164948D-02 -0. 928933D -02 

Q -0. 220778D+02 0.543324D-02 -0. 135074D+02 

THETA 0. OOOOOOD+OO O.OOOOOOD+OO O.OOOOOOD+OO 

VEL -0. 105186D+02 0. 34281 7D+02 -0. 155832D+02 


MATRIX D USING THE FORMULATION OF THE STATE EQUATION 
DX/DT = A*X+B*U+D*V 


FOR CASE # 1 

LINEARIZER DEMONSTRATION AND TEST CASES 


FOR THE PROJECT: 
USER'S GUIDE 


DERIVATIVES 
WRT 
TIME OF: 

X-BODY AXIS FORCE 

Y-BODY AXIS FORCE 

Z-BODY AXIS FORCE 

PITCHING MOMENT 

ROLLING MOMENT 

ALPHA 

-0. 343642D-07 

O.OOOOOOD+OO 

0. 737378D -06 

O.OOOOOOD+OO 

O.OOOOOOD+OO 

Q 

0.1131 92D-06 

O.OOOOOOD+OO 

-0. 242885D -05 

0.605694D -05 

O.OOOOOOD+OO 

THETA 

O.OOOOOOD+OO 

O.OOOOOOD+OO 

O.OOOOOOD+OO 

O.OOOOOOD+OO 

O.OOOOOOD+OO 

VEL 

0. 714203D -03 

0. 3984 92D -06 

0. 332842D-04 

O.OOOOOOD+OO 

O.OOOOOOD+OO 


DERIVATIVES YAWING MOMENT 

WRT 
TIME OF: 

ALPHA O.OOOOOOD+OO 

Q O.OOOOOOD+OO 

THETA O.OOOOOOD+OO 

VEL O.OOOOOOD+OO 


MATRIX H USING THE FORMULATION OF THE OBSERVATION EQUATION: 

Y = H * X + F * U + E * V 

FOR CASE # 1 

LINEARIZER TEST AND DEMONSTRATION CASES 

FOR THE PROJECT: 

USER'S GUIDE 

ALPHA Q THETA VEL 

AN 0. 351 752D+02 -0. 191245D-05 0.150046D-02 0.640771D-02 

AY O.OOOOOOD+OO O.OOOOOOD+OO -0. 150534D -01 0 . 2742 4 80 -09 


MATRIX F USING THE FORMULATION OF THE OBSERVATION EQUATION: 
y = H * X + F * U + E * V 

FOR CASE # 1 

LINEARIZER TEST AND DEMONSTRATION CASES 

FOR THE PROJECT: 

USER'S GUIDE 

ELEVATOR THROTTLE SPEED BRAKE 

AN 0. 41284 5D +01 -0. 180978D -02 0.291699D+00 

AY O.OOOOOOD+OO O.OOOOOOD+OO O.OOOOOOD+OO 
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MATRIX E USING THE FORMULATION OF THE OBSERVATION EQUATION: 
Y=H*X+F*U+E*V 


FOR CASE # 1 

LINEAR IZER DEMONSTRATION AND TEST CASES 

FOR THE PROJECT: 

USER'S GUIDE 

X-BODY AXIS FORCE Y-BODY AXIS FORCE Z-BODY AXIS FORCE PITCHING MOMENT 

AN -0, 377037D -07 O.OOOOOOD+OO -0.214132D-04 O.OOOOOOD+OO 

AY O.OOOOOOD+OO 0.222222D-04 O.OOOOOOD+OO O.OOOOOOD+OO 


YAWING MOMENT 


AN O.OOOOOOD+OO 

AY O.OOOOOOD+OO 


TRIM CONDITIONS FOR CASE # 2 

LINEARIZER TEST AND DEMONSTRATION CASES 

FOR THE PROJECT: 

USER'S GUIDE 

STRAIGHT AND LEVEL TRIM WHILE VARYING ALPHA 
TRIM ACHIEVED 


COEFFICIENT OF LIFT 


= 

0.13221 

COEFFICIENT OF DRAG 


= 

0.00895 

LIFT 

(LBS) 

= 

44376.86258 

DRAG 

(LBS) 

= 

3004.93778 

ALTITUDE 

(FT) 

= 

20000. 

MACH 


= 

0.90000 

VELOCITY 

(FT/SEC) 

= 

933.23196 

EQUIVALENT AIRSPEED 

(KTS ) 

= 

403.42303 

SPEED OF SOUND 

(FT/SEC) 

= 

1036.92440 

GRAVITATIONAL ACCEL 

(FT/SEC**2 ) 

= 

32. 11294 

NORMAL ACCELERATION 

(G-S) 

= 

0.98523 

LOAD FACTOR 

= 

0.98803 

DYNAMIC PRESSURE 

(LB/FT**2 ) 

= 

552.05302 

DENSITY 

(SLUG/FT**3) 

= 

0.00126774 

WEIGHT (^ALTITUDE) 

(LBS) 

= 

44914.60434 

BETA 

(DEG) 

= 

0.00000 

ALPHA 

(DEG) 

= 

-0. 72565 

PHI 

(DEG) 

- 

0.00000 

THETA 

(DEG) 

= 

9.27435 

ALTITUDE RATE 

(FT/SEC) 

= 

162.05403 

GAMMA 

(DEG) 

= 

10.00000 

ROLL RATE 

(DEG/SEC) 

K 

0.00000 

PITCH RATE 

(DEG/SEC) 

= 

0.00000 

YAW RATE 

(DEG/SEC) 


0.00000 

THRUST 

(LBS) 


10804.39673 

SUM OF THE SQUARES 


= 

0.00000 

TRIM PARAMETERS 

TRIM PITCH AXIS PARAMETER 

S 

-0.79364 

TRIM ROLL AXIS PARAMETER 

= 

0.00000 

TRIM YAW AXIS PARAMETER 

= 

0.00000 

TRIM THRUST PARAMETER 


s 

0.22509 


ROLLING MOMENT 

O.OOOOOOD+OO 

O.OOOOOOD+OO 


86 


CONTROL VARIABLES 


ELEVATOR 
THROTTLE 
SPEED BRAKE 


0.06377 

0.22509 

0.00000 


OBSERVATION VARIABLES 


AN 

AY 


0.98522771 GS 
0.00000000 GS 


NON-DIMENSIONAL STABILITY AND CONTROL DERIVATIVES FOR CASE #2 
LINEARIZER TEST AND DEMONSTRATION CASES 


FOR THE PROJECT: 
USER'S GUIDE 




ROLLING 

PITCHING 

YAWING 



SIDE 



MOMENT 

MOMENT 

MOMENT 

DRAG 

LIFT 

FORCE 

ZERO COEFFICIENTS 

-9.37149D-20 

4. 22040D -02 

-6*115950 -1 9 

1.08760D-02 

1 . 57360D -01 

4. 31432D-19 

ROLL RATE 

(RAD/SEC) 

-2. OOOOOD -01 

0. OOOOOD+OO 

-3. 37210D -02 

0. OOOOOD+OO 

0. OOOOOD+OO 

0. OOOOOD+OO 

PITCH RATE 

(RAD/SEC) 

0. OOOOOD+OO 

3. 89530D+00 

0. OOOOOD+OO 

0. OOOOOD+OO 

-1 . 72320D+01 

0. OOOOOD+OO 

YAW RATE 

(RAD/SEC) 

1. 50990D -01 

0. OOOOOD+OO 

-4.04710D -01 

0. OOOOOD+OO 

0. OOOOOD+OO 

0. OOOOOD+OO 

VELOCITY 

(FT/SEC) 

8. 281 1 ID -17 

-1. 01777D-09 

7. 42372D -16 

0. OOOOOD+OO 

7. 18527D -09 

0. OOOOOD+OO 

MACH NUMBER 


8. 58688D -14 

-1 .05535D -06 

7.69784D -13 

0. OOOOOD+OO 

7.45058D -06 

0. OOOOOD+OO 

ALPHA 

(RAD) 

0. OOOOOD+OO 

-1 . 6882 OD -01 

0. OOOOOD+OO 

3. 72570D -01 

4. 87060D+00 

0. OOOOOD+OO 

BETA 

(RAD) 

-1.33450D-01 

0. OOOOOD+OO 

1. 29960D -01 

0. OOOOOD+OO 

0. OOOOOD+OO 

-9. 74030D -01 

ALTITUDE 

(FT) 

0. OOOOOD+OO 

0. OOOOOD+OO 

0. OOOOOD+OO 

0. OOOOOD+OO 

0. OOOOOD+OO 

0. OOOOOD+OO 

ALPHA DOT 

(RAD/SEC) 

0. OOOOOD+OO 

-1 • 18870D+O1 

0. OOOOOD+OO 

0. OOOOOD+OO 

1 . 72320D+01 

0. OOOOOD+OO 

BETADOT 

(RAD/SEC) 

0. OOOOOD+OO 

0. OOOOOD+OO 

0. OOOOOD+OO 

0. OOOOOD+OO 

0. OOOOOD+OO 

0. OOOOOD+OO 

ELEVATOR 


0. OOOOOD+OO 

-6. 95280D-01 

0. OOOOOD+OO 

4. 38310D -02 

5. 72960D -01 

0. OOOOOD+OO 

THROTTLE 


0. OOOOOD+OO 

0. OOOOOD+OO 

0. OOOOOD+OO 

0. OOOOOD+OO 

0. OOOOOD+OO 

0. OOOOOD+OO 

SPEED BRAKE 


0. OOOOOD+OO 

-4. 17500D-01 

0. OOOOOD+OO 

6.49350D-02 

3. 74920D -02 

0. OOOOOD+OO 


VEHICLE STATIC MARGIN IS 3.5% MEAN AERODYNAMIC CHORD STABLE 

AT THIS FLIGHT CONDITION. 


MATRIX A USING THE FORMULATION OF THE STATE EQUATION: 
DX/DT = A * X + B * U + D * V 
FOR CASE # 2 

LINEARIZER TEST AND DEMONSTRATION CASES 


FOR THE PROJECT: 
USER'S GUIDE 


TIME 

DERIVATIVES 

WRT: 

ALPHA 

Q 

THETA 

VEL 

ALPHA 

-0. 12 0900D+01 

0. 100000D+01 

-0.575730D-02 

-0. 701975D -04 

Q 

-0. 1491 89D+01 

-0.221451D+01 

0. 189640D-01 

0. 231368D -03 

THETA 

0. OOOOOOD+OO 

0. lOOOOOD+Ol 

0. OOOOOOD+OO 

0. OOOOOOD+OO 

VEL 

-0. 57686 8D +02 

0. OOOOOOD+OO 

-0. 31 6251D+02 

-0. 46043 5D -02 


MATRIX B USING THE FORMULATION OF THE STATE EQUATION: 
DX/DT = A * X + B * U + D * V 
FOR CASE # 2 

LINEARIZER TEST AND DEMONSTRATION CASES 

FOR THE PROJECT: 

USER'S GUIDE 
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TIME 

DERIVATIVES 
WRT : 

ALPHA 

Q 

THETA 

VEL 


ELEVATOR 


-0. 141 961D+00 
-0.220778D+O2 
0. OOOOOOD+OO 
-0, 105186D+02 


THROTTLE 


0* 448742D -03 
-0. 14781 2D -02 
O.OOOOOOD+OO 
0. 3431 62D+02 


SPEED BRAKE 


-0 . 92 8932D -02 
-0.135074D+02 
0 * 0000000+00 
-0. 155832D+02 


MATRIX D USING THE FORMULATION OF THE STATE EQUATION: 

DX/DT=A*X+B*U+D*V 
FOR CASE #2 

LINEARIZER DEMONSTRATION AND TEST CASES 

FOR THE PROJECT: 

USER'S GUIDE 

X-BODY AXIS FORCE Y-BODY AXIS FORCE Z-BODY AXIS FORCE PITCHING MOMENT ROLLING MOMENT 

DERIVATIVES 

WRT 

TIME OF: 


ALPHA 

0. 934880D -08 

O.OOOOOOD+OO 

0.738119D-06 

O.OOOOOOD+OO 

O.OOOOOOD+OO 

Q 

-0. 307941D-07 

O.OOOOOOD+OO 

-0. 243129D-05 

0. 605694D -05 

O.OOOOOOD+OO 

THETA 

O.OOOOOOD+OO 

O.OOOOOOD+OO 

O.OOOOOOD+OO 

O.OOOOOOD+OO 

O.OOOOOOD+OO 

VEL 

0. 71492 OD -03 

0.49582 9D-13 

-0. 905497D-05 

O.OOOOOOD+OO 

O.OOOOOOD+OO 

DERIVATIVES 

YAWING MOMENT 





WRT 

TIME OF: 






ALPHA 

Q.OGOOGOD+OO 





Q 

O.OOOOOOD+OO 





THETA 

O.OOOOOOD-KIO 





VEL 

O.OOOOOOD+OO 






MATRIX H USING THE 

FORMULATION OF 

THE OBSERVATION I 

EQUATION: 



Y 

= H * X + F * 

■ U + E * V 




FOR CASE # 2 






LINEARIZER TEST AND DEMONSTRATION 

CASES 




FOR THE PROJECT: 
USER'S GUIDE 






ALPHA 

Q 

THETA 

VEL 


AN 

0. 350424D+02 

-0. 128333D -06 

-0.632324D-02 

0. 203434D-02 


AY 

O.OOOOOOD+OO 

O.OOOOOOD+OO 

O.OOOOOOD+OO 

O.OOOOOOD+OO 



MATRIX F USING THE 

FORMULATION OF 

THE OBSERVATION EQUATION: 



Y 

= H * X + F * 

r U + E * V 




FOR CASE # 2 






LINEARIZER TEST AND DEMONSTRATION CASES 




FOR THE PROJECT: 
USER'S GUIDE 






ELEVATOR 

THROTTLE 

SPEED BRAKE 



AN 

0.411323D+01 

0.492845D-03 

0.263288D+00 



AY 

O.OOOOOOD+OO 

O.OOOOOOD+OO 

O.OOOOOOD+OO 
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MATRIX E USING THE FORMULATION OF THE OBSERVATION EQUATION: 
Y=H*X+F*U+E*V 


FOR CASE # 2 

LINEARIZER DEMONSTRATION AND TEST CASES 

FOR THE PROJECT: 

USER'S GUIDE 

X-BODY AXIS FORCE Y-BODY AXIS FORCE Z-BODY AXIS FORCE PITCHING MOMENT ROLLING MOMENT 
AN 0. 102676D -07 O.OOOOOOD+OO -0. 214116D-04 O.OOOOOOD+OO O.OOOOOOD+OO 

AY O.OOOOOOD+OO 0.222222D-04 O.OOOOOOD+OO O.OOOOOOD+OO O.OOOOOOD+OO 


YAWING MOMENT 

AN O.OOOOOOD+OO 

AY O.OOOOOOD+OO 


Example printer output file 2 (unit 2) 


TRIM CONDITIONS FOR CASE # 1 
LINEARIZER TEST AND DEMONSTRATION CASES 

FOR THE PROJECT: 

USER'S GUIDE 

LEVEL TURN WHILE VARYING ALPHA 
TRIM ACHIEVED 


COEFFICIENT OF LIFT 


S 

0.40144 

COEFFICIENT OF DRAG 


a 

0.03058 

LIFT 

(LBS) 

a 

134741.68807 

DRAG 

(LBS) 

* 

10265.70661 

ALTITUDE 

(FT) 

s 

20000. 

MACH 


a 

0.90000 

VELOCITY 

(FT/ SEC) 

= 

933.23196 

EQUIVALENT AIRSPEED 

(KTS) 

= 

403.42303 

SPEED OF SOUND 

(FT/SEC) 

= 

1036.92440 

GRAVITATIONAL ACCEL 

(FT/SEC**2 ) 

a 

32.11294 

NORMAL ACCELERATION 

(G-S) 

a 

3.00163 

LOAD FACTOR 

a 

2.99995 

DYNAMIC PRESSURE 

(LBS/FT**2 ) 

a 

552.05302 

DENSITY 

(SLUG/FT**3) 

s 

0.00126774 

WEIGHT (^ALTITUDE) 

(LBS) 

= 

44914.60434 

BETA 

(DEG) 

= 

0.03193 

ALPHA 

(DEG) 

a 

2.66824 

PHI 

(DEG) 

= 

70.62122 

THETA 

(DEG) 

a 

0.91607 

ALTITUDE RATE 

(FT/ SEC) 

= 

0.00000 

GAMMA 

(DEG) 

s 

0.00000 

ROLL RATE 

(DEG/SEC) 

a 

-0.08951 

PITCH RATE 

(DEG/SEC) 

= 

5.28086 

YAW RATE 

(DEG/SEC) 

= 

1.85749 

THRUST 

(LBS) 

= 

10277.03515 

SUM OF THE SQUARES 


a 

0.00000 

TRIM PARAMETERS 

TRIM PITCH AXIS PARAMETER 

a 

-0.66958 

TRIM ROLL AXIS PARAMETER 

= 

-0.01526 

TRIM YAW AXIS PARAMETER 

s 

-0.02125 

TRIM THRUST PARAMETER 


= 

0.21410 

CONTROL VARIABLES 

ELEVATOR 


a 

0.05380 

THROTTLE 


a 

0.21410 

SPEED BRAKE 

(DEG) 

= 

0.00000 
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OBSERVATION VARIABLES 


3.00163339 GS 
0 . 94136286 GS 

VEHICLE STATIC MARGIN IS 3.5% MEAN AERODYNAMIC CHORD STABLE 
AT THIS FLIGHT CONDITION 


AN 

AY 


TRIM CONDITIONS FOR CASE # 2 

LINEARIZER TEST AND DEMONSTRATION CASES 

FOR THE PROJECT: 

USER'S GUIDE 


STRAIGHT AND LEVEL TRIM WHILE VARYING ALPHA 
TRIM ACHIEVED 


COEFFICIENT OF LIFT 


= 

0.13221 

COEFFICIENT OF DRAG 


ss 

0.00895 

LIFT 

(LBS) 

= 

44376.86258 

DRAG 

(LBS) 

= 

3004.93778 

ALTITUDE 

(FT) 

= 

20000. 

MACH 


= 

0.90000 

VELOCITY 

(FT/ SEC) 

= 

933.23196 

EQUIVALENT AIRSPEED 

(KTS) 

= 

403.42303 

SPEED OF SOUND 

(FT/SEC) 

* 

1036.92440 

GRAVITATIONAL ACCEL 

(FT/SEC**2 ) 

= 

32.11294 

NORMAL ACCELERATION 

(G-S) 

= 

0.98523 

LOAD FACTOR 


= 

0.98803 

DYNAMIC PRESSURE 

(L8/FT**2 ) 

= 

552.05302 

DENSITY 

(SLUG/FT**3) 

= 

0.00126774 

WEIGHT (^ALTITUDE) 

(LBS) 

= 

44914.60434 

BETA 

(DEG) 

= 

0.00000 

ALPHA 

(DEG) 

= 

-0.72565 

PHI 

(DEG) 

- 

0.00000 

THETA 

(DEG) 

= 

9.27435 

ALTITUDE RATE 

(FT/ SEC) 

= 

162.05403 

GAMMA 

(DEG) 

= 

10.00000 

ROLL RATE 

(DEG/SEC) 

= 

0.00000 

PITCH RATE 

(DEG/SEC) 

= 

0.00000 

YAW RATE 

(DEG/SEC) 

= 

0.00000 

THRUST 

(LBS) 

= 

10804.39673 

SUM OF THE SQUARES 



0.00000 

TRIM PARAMETERS 




TRIM PITCH AXIS PARAMETER 

- 

-0.79364 

TRIM ROLL AXIS PARAMETER 

= 

0.00000 

TRIM YAW AXIS PARAMETER 

= 

0.00000 

TRIM THRUST PARAMETER 


= 

0.22509 

CONTROL VARIABLES 




ELEVATOR 


— 

0.06377 

THROTTLE 


= 

0.22509 

SPEED BRAKE 


= 

0.00000 

OBSERVATION VARIABLES 




AN 


= 

0.98522771 

AY 


= 

0.00000000 


VEHICLE STATIC MARGIN IS 3.5% MEAN AERODYNAMIC CHORD STABLE 
AT THIS FLIGHT CONDITION 


APPENDIX Is EXAMPLE USER-SUPPLIED SUBROUTINES 


The following subroutines are examples of user-supplied routines that provide 
the aerodynamic, control, and engine models to LINEAR. These subroutines are based 
on an F-1 5 aircraft simulation and are typical of the routines needed to interface 
LINEAR to a set of nonlinear simulation models. These subroutines are meant to 
illustrate the use of the named common blocks to communicate between LINEAR and 
the user’s routines. These subroutines are used with all examples in this report. 
Included with this report are microfiche listings of these subroutines. 


Aerodynamic Model Subroutines 

The following two subroutines define a linear aerodynamic model. Even though 
this model is greatly simplified from the typical nonlinear aerodynamic model, the 
example illustrates the functions of the subroutines ADATIN and CCALC. 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 


C 

C 

C 

C 

C 

C 


C 

C 

C 

c 


c 

c 

c 

c 

c 

c 

c 


c 


SUBROUTINE ADATIN 

.EXAMPLE AERODYNAMIC DATA DEFINITION (OR INPUT) SUBROUTINE 

•ROUTINE TO DEFINE STABILITY AND CONTROL DERIVATIVES FOR THE 
AERODYNAMIC MODEL 


COMMON BLOCK TO COMMUNICATE AERODYNAMIC DATA BETWEEN THE 
SUBROUTINES ADATIN AND CCALC 


CLB 

,CLP 

,CLR 

> 

CLDA 

,CLDR 

,CLDT 

> 

CM0 

,CMA 

,CMQ 

,CMAD 

CMDE 

,CM$B 

» 


CNB 

,CNP 

,CNR 

> 

CNDA 

,CNDR 

,CNDT 

> 

CD0 

,CDA 

,CDDE 

,CDSB 

CLFT0 

.CLFTA 

,CLFTQ 

,CLFTAD 

CLFTDE 

,CLFTSB 

» 


CYB 

,CYDA 

,CYDR 

,CYDT 


.ROLLING MOMENT DERIVATIVES. 

STABILITY DERIVATIVES WITH RESPECT TO 
SIDESLIP, ROLL RATE, YAW RATE . . . 

CLB = -1.3345 E-01 
CLP = -2.0000 E-01 
CLR = 1.5099 E-01 

CONTROL DERIVATIVES WITH RESPECT TO 
AILERON, RUDDER, DIFFERENTIAL TAIL . . . 

CLDA = 2.6356 E-02 
CLDR = -2.3859 E-03 
CLDT = 4.0107 E-02 


.PITCHING MOMENT DERIVATIVES. 


STABILITY DERIVATIVES WITH RESPECT TO ANGLE OF 
ATTACK, PITCH RATE, ANGLE OF ATTACK RATE . . . 


CM0 

= 4.2204 

E-02 

CMA 

= -1.6882 

E-01 

CMQ 

= 3.8953 


CMAD 

= -1.1887 

E+01 
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CONTROL DERIVATIVES WITH RESPECT TO 
ELEVATOR, SPEED BRAKE . . . 

CMDE = -6.9528 E-01 
CMSB = -4.1750 E-01 

....YAWING MOMENT DERIVATIVES. 

STABILITY DERIVATIVES WITH RESPECT TO 
SIDESLIP, ROLL RATE, YAW RATE . . . 

CNB = 1.2996 E-01 

CNP = -3.3721 E-02 
CNR = -4.0471 E-01 

CONTROL DERIVATIVES WITH RESPECT TO 
AILERON, RUDDER, DIFFERENTIAL TAIL . . . 

CNDA = 2.1917 E-03 
CNDR = -6.9763 E-02 
CNDT = 3.0531 E-02 

....COEFFICIENT OF DRAG DERIVATIVES. 

STABILITY DERIVATIVES WITH RESPECT TO 
ANGLE OF ATTACK . . . 

CDO = 1.0876 E-02 

CDA = 3.7257 E-01 


CONTROL DERIVATIVES WITH RESPECT TO 
ELEVATOR, SPEED BRAKE . . . 

CDDE = 4.3831 E-02 
CDSB = 6.4935 E-02 

....COEFFICIENT OF LIFT DERIVATIVES. 

STABILITY DERIVATIVES WITH RESPECT TO ANGLE OF 
ATTACK, PITCH RATE, ANGLE OF ATTACK RATE . . . 

CLFTO = 1.5736 E-01 

CLFTA = 4.8706 
CLFTQ = -1.7232 E+Ol 
CLFTAD= 1.7232 E+Ol 

CONTROL DERIVATIVES WITH RESPECT TO 
ELEVATOR, SPEED BRAKE . . . 

CLFTDE= 5.7296 E-01 
CLFTSB= 3.7492 E-02 

SIDEFORCE COEFFICIENT DERIVATIVES. 

STABILITY DERIVATIVES WITH RESPECT TO 
SIDESLIP . . . 

CYB = -9.7403 E-01 

CONTROL DERIVATIVES WITH RESPECT TO 
AILERON, RUDDER, DIFFENTIAL TAIL 

CYDA = -1.1516 E-03 
CYDR * -1.5041 E-01 
CYDT = -7.9315 E-02 
C 

RETURN 

END 
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SUBROUTINE CCALC 


C 

c.. 

c 

c.. 

c 

c 

c 

c 


c 


c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 


c 


c 

c. 

c 


c 

c. 

c 


c. 

c 


c 

c. 


EXAMPLE AERODYNAMIC MODEL, 

ROUTINE TO CALCULATE THE AERODYNAMIC FORCE AND MOMENT COEFFICIENTS. 

COMMON BLOCKS CONTAINING STATE, CONTROL, AND AIR 
DATA PARAMETERS 

COMMON /DR VO UT/ F(13),DF(13) 

COMMON /CONTRL/ DC(30) 

COMMON /DATAIN/ S,B,CBAR,AMSS,AIX,AIY,AIZ,AIXZ ,AIXY,AIYZ,AIXE 
COMMON /TRIGFN/ SINALP,COSALP,SINBTA,COSBTA,SINPH I , COS PHI , 

S I N PS I , COS PS I , S I NTHA , COSTHA 
COMMON /SIMOUT/ AMCH,QBAR,GMA,DEL,UB, VB,WB, VEAS, VCAS 
COMMON /CGSHFT/ DELX,DELY,DELZ 

COMMON BLOCK TO OUTPUT AERODYNAMIC FORCE AND MOMENT 
COEFFICIENTS 

COMMON /CLCOUT/ CL ,CM ,CN ,CD ,CLFT ,CY 

COMMON BLOCK TO COMMUNICATE AERODYNAMIC DATA BETWEEN 
THE SUBROUTINES ADATIN AND CCALC 

COMMON /ARODAT/ CLB ,CLP ,CLR , 

CLDA ,CLDR ,CLDT , 

CMO ,CMA ,CMQ ,CMAD , 

CMDE ,CMSB , 

CNB ,CNP ,CNR , 

CNDA ,CNDR ,CNDT , 

CDO ,CDA ,CDDE ,CDSB , 

CLFTO ,CLFTA ,CLFTQ ,CLFTAD, 

CLFTDE.CLFTSB, 

CYB ,CYDA ,CYDR ,CYDT 

EQUIVALENCE VARIABLE NAMES 

EQUIVALENCE (T , F( 1)), 

(P , F( 2)),(Q , F( 3) ) ,(R , F( 4)), 

(V , F( 5)), (ALP , F( 6 ) ) , (BTA , F( 7)), 

(THA , F( 8) ) ,(PS I , F( 9)), (PHI , F(10)), 

(H , F(ll ) ) ,(X , F( 12 ) ) , ( Y , F ( 1 3 ) ) , 

(TDOT ,DF( 1)), 

(PDOT t DF ( 2 ) ) , (QDOT ,DF( 3)),(RD0T ,DF( 4)), 

(VDOT , DF ( 5) ) ,(ALPDOT,DF( 6) ) , (BTADOT ,DF ( 7)), 
(THADOT ,DF ( 8) ) , (PS IDOT,DF( 9 ) ) , (PHIDOT.DF(IO) ) , 
(HDOT ,DF (11 ) ) , (XDOT ,DF( 12 ) ) ,(YDOT ,DF(13)) 

EQUIVALENCE (DA ,DC( 1)),(DE ,DC( 5)),(DT ,DC( 8)), 

(DR ,DC( 9 ) ) ,(DSB ,DC(10)) 

.COMPUTE TERMS NEEDED WITH ROTATIONAL DERIVATIVES. 


V2 

= 2.0 

*V 

B2V 

=B 

/V 2 

C2V 

=CBAR 

/V 2 


.ROLLING MOMENT COEFFICIENT. 

CL =CLB * BTA 4CLDA *DA 4CLDR *DR 4CLDT *DT 

462V *(CLP *P 4CLR *R ) 

.PITCHING MOMENT COEFFICIENT. 

CM <MO + CMA *ALP 4CMDE *DE 4CMSB *DSB 

. 4C2V *(CMQ *Q 4CMAD *ALPDOT ) 

.YAWING MOMENT COEFFICIENT. 
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CN <NB 

* BTA 

4CNDA 

*DA 

■CNDR 

*DR 

-k:ndt 

*DT 

r 

+B2V 

*(CNP 

*P 

-k:nr 

*R 

) 



L 

c. 

r 

...COEFFICIENT 

OF DRAG. 







r 

CD =CD0 

+ CDA 

*ALP 

-ODE 

*DE 

+CDSB 

*DSB 


c. 

r 

...COEFFICIENT 

OF LIFT. 







L 

CLFT =CLFT0 

+ CLFTA 

*ALP 

-OFTDE*DE 

+CLFTSB*DSB 


r 

-ov 

*(CLFTQ 

*Q 

-OFTAD*ALPDOT 

) 



\j 

c. 

r 

...SIDEFORCE COEFFICIENT. 







L 

CY =C YB 

*BTA 

4CYDA 

*DA 

+C YDR 

*DR 

-k:ydt 

*DT 


RETURN 

END 


Engine Model Interface Subroutine 

The following subroutine, IFENGN, both defines an engine model and provides the 
interface to LINEAR. in normal usage this subroutine would call subroutines that 
provide the detailed thrust, engine rotation, and fuel consumption modeling; the 
information from these subroutines would be transferred into the /ENGSTF/ common 
block. 


SUBROUTINE IFENGN 
C 

C EXAMPLE SUBROUTINE TO PROVIDE PROPULSION SYSTEM MODEL. 

C 

C ROUTINE TO COMPUTE PROPULSION SYSTEM INFORMATION FOR LINEAR 

C 

C THIS SUBROUTINE IS THE INTERFACE BETWEEN THE DETAILED ENGINE 
C MODELING SUBROUTINES AND LINEAR. 

C 

C INPUT COMMON BLOCK CONTAINING INFORMATION ON THRUST 

C REQUEST TO ENGINE 

C 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 


COMMON /CONTRL/ DC (30) 

OUTPUT COMMON BLOCK CONTAINING DETAILED INFORMATION 
ON EACH OF UP TO FOUR SEPARATE ENGINES 

..ROUTINE TO COMPUTE ENGINE PARAMETERS: 


THRUST (I ) 
TLOCAT (1,0) 
XYANGL (I ) 

XZANGL (I ) 

TVANXY (I ) 

TVANXZ(I) 

DXTHRS (I ) 

EIX(I) 


THRUST CREATED BY EACH ENGINE 
LOCATION OF EACH ENGINE IN THE X-Y-Z PLANE 
ANGLE IN X-Y BODY AXIS PLANE AT WHICH EACH 
ENGINE IS MOUNTED 

ANGLE IN X-Z BODY AXIS PLANE AT WHICH EACH 
ENGINE IS MOUNTED 

ANGLE IN THE X-Y ENGINE AXIS PLANE OF THE 
THRUST VECTOR 

ANGLE IN THE X-Z ENGINE AXIS PLANE OF THE 
THRUST VECTOR 

DISTANCE BETWEEN THE ENGINE C.G. AND THE 
THRUST POINT 

ROTATIONAL INERTIA OF EACH ENGINE 


C AMSENG(I) = MASS OF EACH ENGINE 

C ENGOMG(I) = ROTATIONAL VELOCITY (RAD/SEC) 

C 

COMMON /ENGSTF/ THRUST (4 ) ,TLOCAT (4, 3),XYANGL(4) ,XZANGL(4) , 
TVANXY (4 ),TVANXZ(4) ,DXTHRS (4 ) , 

. E I X (4 ) ,AMSENG( 4),ENG0MG(4) 

C 

C EQUIVALENCE VARIABLE NAMES 

C 

EQUIVALENCE (THR ,DC(12)) 

C 

C ASSUME THRUST PER ENGINE IS HALF VEHICLE WEIGHT. 

C 

THR UST(1 )= 24000.0 *THR 
THRUST(2)= 24000.0 *THR 
C 

C LET ALL OTHER PARAMETERS DEFAULT TO ZERO 

C 

RETURN 

END 


Control Model Subroutine 

The subroutine UCNTRL provides an interface between the trim parameters and the 
surface deflections • Figure 8 illustrates the gearing model implemented in the 
example UCNTRL subroutine. The thrust demand parameter is also set and passed to 
the subroutine IFENGN. 



UCNTRL 


Figure 8 . Gearing model in example 
UCNTRL subroutine . 
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SUBROUTINE UCNTRL 
C 

C. .. .EXAMPLE TRIM/CONTROL SURFACE INTERFACE ROUTINE. 

C 

C 

C.... ROUTINE TO CONVERT TRIM INPUTS INTO CONTROL SURFACE DEFLECTIONS. 

C 

C INPUT COMMON BLOCK CONTAINING TRIM PARAMETERS 

C 

COMMON /CTPARM/ DES ,DAS ,DRS .THRSTX 
C 

C OUTPUT COMMON BLOCK CONTAINING CONTROL SURFACE DEFLECTIONS 

C 

COMMON /CONTRL/ DC (30) 

C 

C EQUIVALENCE VARIABLE NAMES. 

C 

EQUIVALENCE (DA ,DC( 1)),(DE ,DC( 5)),(DT ,DC( 8)), 

(DR ,DC( 9 ) ) ,(DCB ,DC(10) ) ,(THR ,DC(12)) 

C 

DATA DGR /57. 29578/ 

C 

C.... CONVERT FROM INCHES OF STICK AND PEDAL TO DEGREES OF SURFACE 
C DEFLECTION 
C 

DA =DAS *( 20.0 / 4.0 ) 

DE =DES *(-25.0 / 5.43) 

DR =DRS *( 30.0 / 3.25) 

C 

C.. . .SET DIFFERENTIAL TAIL BASED ON AILERON COMMAND. 

C 

DT =DA / 4.0 

C 

C.... CONVERT THRUST TRIM PARAMETER TO PERCENT THROTTLE COMMAND. 

C 

THR = 0.0 

IF (THRSTX.GE. 0.0 ) THR =THRSTX 
C 

C....USE SPEED BRAKE IF NEEDED. 

C 

DSB =0.0 

IF(THRSTX.LT. 0.0 ) DSB =THRSTX*(-45.0 ) 

C 

C... .CONVERT SURFACE COMMANDS TO RADIANS. 

C 


DA 

=DA 

/DGR 

DE 

=DE 

/DGR 

DR 

=DR 

/DGR 

DT 

=DT 

/DGR 

DSB 

=DSB 

/DGR 


C 

RETURN 

END 


Mass and Geometry Model Subroutine 

The following subroutine, MASGEO, is an example of the mass and geometry sub- 
routine that the user must provide to LINEAR. If the mass and geometry character- 
istics of the aircraft do not change or can be easily specified as a function of 
flight condition and surface setting (vehicle configuration), then the mass and 
geometry characteristics of the vehicle can be defined using the input file, and 
MASGEO would be a dummy subroutine, as is the following example. However, for an 
aircraft in which these mass and geometry properties are complicated or already 
defined using FORTRAN subroutines, MASGEO provides an interface between LINEAR and 
those subroutines . 
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SUBROUTINE MASGEO 


.... SUBROUTINE TO COMPUTE THE MASS AND GEOMETRY PROPERTIES OF THE 
AIRCRAFT. 

COMMON /CONPOS/ DAP, FLAT, DATRIM.DEP.FLON, DETRIM, DRP.FPED.DRTRIM, 
. DSB P , DF P , PL A PL , PLA PR , THS L , THSR 

COMMON /CONTRL/ DC(30) 

COMMON /DATAIN/ S,B,CBAR,AMSS,AIX,AIY,AIZ,AIXZ,AIXY,AIYZ,AIXE 
COMMON /SIMOUT/ AMCH.QBAR.GMA.DEL.UB, VB.WB, VEAS.VCAS 
COMMON /DRVOUT/ F(13),DF(13) 

EQUIVALENCE(T , F( 1)), 

(P , F( 2 ) ) ,(Q , F( 3 ) ) , (R , F( 4)), 

(V , F( 5)), (ALP , F( 6)),(BTA , F( 7)), 

(THA , F( 8) ) , (PS I , F( 9)), (PHI , F(10)), 

(H , F(11)),(X , F( 1 2 ) ) ,(Y , F( 13 ) ) , 

(TDOT ,DF( 1)), 

. (PDOT ,DF ( 2)),(QDOT ,DF( 3)),(RDOT ,DF( 4)), 

. (VDOT ,DF( 5)),(ALPDOT,DF( 6)),(BTAD0T,DF( 7)), 

. (THADOT,DF( 8)),(PSID0T,DF( 9 ) ) , (PHIDOT.DF ( 10 ) ) , 

(HDOT ,DF(11)),(XD0T ,DF(12)),(YD0T , DF ( 1 3 ) ) , 

RETURN 

END 
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APPENDIX J: REVISIONS TO MICORFICHE SUPPLEMENT 


The following listing of subroutine ENGINE incorporates revisions not contained 
in the microfiche supplement included with this report (Program LINEAR and Subrou- 
tines Utilized by LINEAR) • This listing should be used in place of the subroutine 
ENGINE listing on the microfiche supplement. 
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